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Abstract 


Morphologic  and  hydraulic  data  were  gathered  from  60  gravel 
channels  (20/reach)  on  the  Beauty  Creek  Flats  section  of  the 
proglacial  Sunvapta  River  during  August  4-15,  1977-  Four 
power  functions  (log-linear  ;  width,  average  depth,  average 
surface  velocity  and  water  surface  slope)  are  determined  for 
each  reach  and  for  Beauty  Creek  Flats  as  a  whole  with 
discharge  as  the  independent  parameter  in  all  cases-  A 
downstream  hydraulic  geometry  interpretation  indicates  that 
the  downstream  reduction  in  the  median  grain  size,  D^q ,  is 
mainly  responsible  for  the  variation  in  channel  behaviour. 

In  all  three  reaches  of  Beauty  Creek  Flats  the  primary 
adjustment  of  the  channels  to  an  increased  discharge  is  an 
increase  in  the  width  of  the  flow. 

It  appears  that  the  width  exponent  (b)  may  be  used  as 
an  index  of  braiding  intensity. 

The  b-f-m  diagram  of  Rhodes  (1977)  is  used  to  determine 
the  at-a-station  variation  in  channel  behaviour  downstream. 

The  hydraulic  geometry  relations  determined  for  Beauty 
Creek  Flats  as  a  whole  are  compared  to  corresponding 
relations  from  other  gravelly,  glac iof luvial  environments 
(Fahnestock  (1963),  Church  (1972)  and  Arnborg  (1955)).  The 
best  comparison  is  between  the  Sunwapta  River  equations  and 
those  of  Fahnestock  determined  for  the  White  River.  A 
difference  in  channel  behaviour  between  the  Sunwapta  River 
and  the  White  River  as  a  group  and  the  Baffin  Island  sandar 
(Church)  and  the  Hof f ellssandur  (Arnborg)  as  another  group 
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is  suggested  to  be  caused  by  greater  bank  stability  in  the 
Arctic  studies. 

It  is  suggested  that  there  is  a  similarity  in  channel 
response  among  areas  of  braided  gravel  outwash  provided 
associated  climatic  conditions  are  similar. 

The  width  and  velocity  relations  of  Church  and  the 
depth  relation  of  this  study  (eq.  4.14)  are  suggested  for 
estimating  channel  response  on  braided  gravel  outwash  in  the 
Arctic  environment.  The  width  and  velocity  relation  of 
Fahnestock  and  the  depth  relation  of  this  study  are 
suggested  for  similar  usage  in  the  non-Arctic  environment. 
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NOTATION 


The  following  points  are  to  be  noted: 

1.  The  use  of  an  overbar  in  the  text  indicates  a  mean  value. 

2.  Some  symbols  have  more  than  one  meaning.  The  context  of 
its  usage  will  indicate  which  meaning  is  appropriate. 
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MEANING(S) 

DIMENSION (S) 

b  J  f  y  Hi  )  z 

exponents  of  the  hydraulic  geometry 
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coefficients  of  the  hydraulic 
geometry  relations  for  width,  depth, 
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specific  weight  =  Pg  ;  Yf  _  specific 
weight  of  fluid;  -  specific  weight 
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fluid  shear  stress;  tq  -  average 
bottom  shear  stress;  x^  -  channel 
bank  shear  stress;  xc  -  critical 
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M  L~2 
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theoretical  discharge  correction 
factor 
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dynamic  viscosity  of  water 
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kinematic  viscosity  of  water  =  y/p 
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p 

density;  p^  -  fluid  density  ; 
p  -  sediment  density 
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0.6  of  depth  of  flow*;  V  -  surface 
velocity;  V  -_theoretical  surface 
velocity;  V  -  average  velocity 

of  flow  cross  section 

river  stage;  channel  gradient 

slope  of  the  water  surface  ;  slope 
of  the  channel  bed 
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50  and  90  percentile  of  a  grain 
size  distribution  by  weight 

the  coefficient  of  determination  in 
a  least  squares  regression 
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Mannings  n,  a  total  roughness 
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gravitational  constant 

channel  width  at  the  surface  of  the 
flow 

dimensionless  channel  width  at  the 
surface 

channel  depth 

dimensionless  channel  depth 

corrected  channel  depth 

width  to  corrected  depth  ratio 

hydraulic  radius  of  channel  ; 
apparent  density  of  the  sediment 
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dimensionless  water  discharge 
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1.  Introduction 


1_t_l  The  Topic 

This  investigation  centers  on  the  hydraulic  geometry  of 
natural  river  channels.  Hydraulic  geometry  refers  to  the 
manner  in  which  river  channel  width,  depth,  and  velocity 
vary  in  response  to  a  change  in  discharge  of  the 
sediment- water  complex.  These  three  parameters  can  be 
considered  as  summary  parameters  reflecting  the  net  result 
of  complex  interactions  between  the  numerous  morphologic  and 
hydraulic  factors  involved  in  open  channel  flow  in  natural 
river  channels.  While  hydraulic  geometry  deals  specifically 
with  only  the  relationships  between  width,  depth  and 
velocity  to  discharge  all  other  morphologic  and  hydraulic 
parameters  involved  in  open  channel  flow  are  being 
implicitly  considered  (Rhodes,  1977,  p.  73). 

The  conventional  mathematical  model  for  these 
relationships  is  a  power  function  of  the  form 

b 

Y  =  aX 

The  exponent  (b)  being  the  slope  of  the  relationship 
indicates  the  sense  (  positive  or  negative  correlation  )  and 
rate  of  response  of  the  dependent  parameter  (Y)  to  changes 
in  the  independent  parameter  (X).  The  rate  of  response  of 
the  width,  depth  or  velocity  (Y)  to  an  imposed  discharge  (X) 
is  a  function  of  the  nature  of  the  channel  bed  and  bank 
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materials  which  the  determine  the  shape  of  the  channel  cross 
section.  The  hydraulic  geometry  exponent  is  therefore 
reflecting  how  a  particular  channel  morphology  is  responding 
to  the  hydraulic  characteristics  of  the  flow  it  is 
transmitting. 

Hydraulic  geometry  studies  may  be  separated  into  two 
categories,  namely  at-a-station  hydraulic  geometry  and 
downstream  hydraulic  geometry.  In  the  former  case 
width-discharge,  depth-discharge  and  velocity-discharge 
relations  are  developed  for  a  single  (ideally)  channel  cross 
section  (station)  by  successive  measurements  of  channel 
width,  depth  and  flow  velocity  at  different  river  stages 
corresponding  to  different  discharges.  The  emphasis  is  on 
variations  in  the  shape  of  the  flow  cross  section  as 
indicated  by  the  combined  changes  in  cross  section  width, 
depth  and  velocity.  In  downstream  hydraulic  geometry  studies 
rather  than  combine  the  changes  in  width,  depth  and  velocity 
and  deal  with  the  flow  cross  section  as  a  whole,  changes  in 
these  flow  parameters  are  considered  individually  in  terms 
of  comparing  their  response  to  changing  discharge  at 
successive  downriver  locations  (channel  cross  sections). 

The  discharge  variations  encountered  in  at-a-station 
studies  are  weather  induced  and  may  occur  daily,  weekly  or 
monthly  depending  on  the  frequency  of  the  fluctuations  in 
the  weather  and  on  the  characteristics  of  the  drainage 
basin.  Weather  induced  temporal  changes  in  discharge  are 
also  reflected  in  the  discharges  involved  in  downstream 
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studies  but  the  primary  reason  for  the  change  in  discharge 
in  these  studies  is  due  to  the  discharge  additions  from  the 
major  and  minor  tributaries  joining  a  trunk  channel  at 
irregular  intervals  along  its  length.  As  a  result,  in  most 
rivers  the  total  discharge  will  increase  as  the  distance 
downstream  increases.  In  many  drainage  basins  the  discharge 
of  the  major  channels  are  monitored.  In  these  basins  the 
discharge  at  any  point  along  the  trunk  channel  can  be 
determined  by  summing  the  individual  discharge  contributions 
of  the  upstream  tributaries.  The  correlation  of  changes  in 
width,  depth  and  velocity  with  changes  in  discharge  at  a 
particular  station  on  a  trunk  channel  and  the  comparison  to 
similar  correlations  from  other  stations  along  the  same 
channel  constitutes  downstream  hydraulic  geometry. 

An  important  concept  in  downstream  hydraulic  geometry 
is  that  of  dominant  discharge.  A  dominant  discharge  must  be 
stated  in  any  downstream  hydraulic  geometry  study  since  it 
should  represent  the  discharge  at  which  or  in  the  area  of 
which  the  measurements  of  width,  depth  and  velocity  were 
obtained.  As  such  it  represents  a  type  of  reference  point. 

At  present  there  is  no  consensus  of  opinion  as  to  what  flow 
frequency  should  represent  a  dominant  discharge  (Bray,  1972, 
p.  14).  In  reviewing  the  concept  Bray  (p.  9-14)  presents 
several  definitions  of  dominant  discharge  taken  from  the 
literature;  they  are  as  follows: 
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1.  Inglis  (1S41)  -  the  discharge  which  controls  the  meander 
length  and  breadth;  it  appears  to  be  slightly  less  than 
barkf  ull. 

2.  Leopold  and  Maddock  (1953)  -  the  mean  annual  discharge 
was  used  as  the  representative  discharge  for  the 
determination  of  the  hydraulic  geometry  relations, 
however  they  were  not  implying  that  it  should  be 
considered  as  a  dominant  discharge. 

3.  Nixon  (1959)  -  the  dominant  discharge  was  defined  as  the 

bankfull  discharge- 

4.  Leopold,  Wolman  and  Miller  (1964)  -  they  stated  that  the 
effective  discharge  can  often  be  approximated  by 
bankfull  discharge. 

5.  Prins  (1969)  -  the  dominant  discharge  was  defined  by  the 
amount  of  bedload  transport  it  produced. 

6.  Blench  (1957)  -  the  dominant  discharge  was  defined  as 
that  discharge  which  was  associated  with  the  average 
load  of  a  stream  for  any  long  period  of  time- 

7.  Benson  and  Thomas  (1966)  -  dominant  discharge  was 
defined  as  the  discharge  at  which  the  maximum  suspended 
sediment  load  is  transported. 

Bray  (p.  14)  concludes  his  review  by  stating  that  despite 

the  disparity  over  the  definition  of  a  dominant  discharge  it 
is  generally  accepted  that  a  dominant  discharge  is  a 
formative  discharge,  ie.  a  discharge  which  is  capable  of 
physical  alteration  of  the  channel  cross  section-  In  the 
context  of  the  previous  discussion  where  the  exponent  of  the 
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power  function  was  defined  as  an  indicator  of  the  rate  of 
adjustment  of  the  channel  cross  section  to  an  imposed 
discharge,  the  formative  nature  of  a  dominant  discharge  is 
mandatory. 

Equally  important  to  downstream  hydraulic  geometry 
studies  is  the  concept  of  constant  frequency  of  discharge. 
This  concept  expresses  the  necessity  of  comparing  different 
stations  along  a  river  using  data  that  was  obtained  at 
discharges  of  equal  flow  frequency.  For  instance,  as 
discussed  in  Leopold  and  Maddock's  (1953)  paper,  if  data 
obtained  at  an  upstream  station  was  secured  at  the  mean 
annual  discharge  for  that  station,  then  to  compare  its 
hydraulic  geometry  to  that  of  a  downstream  station  the  data 
on  the  latter  must  have  been  obtained  at  i  mean  annual 
discharge.  The  absolute  value  of  the  mean  annual  discharge 
for  the  downstream  station  will  exceed  that  of  the  upstream 
station  but  the  frequency  of  occurrence  of  the  two  flows 
will  be  the  same. 

Hydraulic  geometry  is  a  tool,  a  useful  means  of 
summarizing  the  very  complicated  interactions  between  the 
flow  and  the  bed  and  bank  materials  of  a  natural  channel 
(Rhodes,  1977,  p.  94).  The  utility  of  hydraulic  geometry 
studies  lies  in  the  predictive  capability  of  the  equations 
which  define  the  response  of  the  channel  and  flow 
parameters.  These  equations  serve  as  indicators  of  channel 
behaviour  as  a  function  of  discharge  and  as  such  become 
important  sources  of  information  to  anyone  concerned  with 


6 


the  response  of  a  particular  river  channel  to  changes  in  its 
discharge. 


1-2  The  Objectives  Of  The  Investigation 

Morphologic  and  hydraulic  data  was  acquired  from  the 
gravel  channels  of  the  study  area  with  the  following 
objectives  in  mind: 

1.  The  characterization  cf  the  morphologic  and  hydraulic 
nature  of  the  active  channels. 

2.  The  determination  of  a  variation  in  channel  behaviour 
through  the  development  of  hydraulic  geometry  relations 
and  the  interpretation  of  the  cause  of  this  variation- 

3.  A  comparison  to  channel  behaviour  in  other  areas  of 
gravelly  braided  outwash. 

4.  The  selection  of  hydraulic  geometry  relations  which 
could  be  used  as  estimators  of  channel  behaviour  in  any 
area  of  gravelly  braided  outwash. 

1 . 3  The  Study  Area 
1.3.1  Location 

The  river  chosen  for  this  study  is  the  Sunwapta, 
located  at  the  southern  end  of  Jasper  National  Park  in 
western  Alberta.  The  river  is  approximately  50  km  in  length, 
joining  with  the  Athabasca  River  at  Athabasca  Falls.  The 
uppermost  20  km  of  the  river  is  discontinuously  braided, 
flowing  through  gravelly  outwash  deposited  proglacially.  The 
lower  9  km  of  the  upper  20  km  constitutes  a  continuous 


stretch  of  braided  gravel  flats  known  locally  as  Beauty 
Creek  Flats.  These  gravel  flats  were  chosen  as  the  study 
area  (figure  1-1,  p.  8). 
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JL=.3-_2  Climate 

The  climate  of  the  region  (Canadian  Rocky  Mountains) 
has  been  described  by  Heusser  (1956)  as  continental  in 
character-  Heusser  (tbl.  2,  p.  268)  presents  the  temperature 
and  precipitation  extremes  from  discontinuous  records 
started  at  Jasper  in  1916.  At  the  time  of  publication  of  his 
paper  the  maximum  and  minimun  recorded  temperatures  were  98 
F°  (36.6  C°)  and  -53  F°  (-47.2  C°)  in  1941  and  1916 
respectively.  The  maximum  and  minimum  precipitation  was 
22.76  in.  (57.86  cm)  and  8.47  in.  (21.51  cm)  in  1948  and 
1941  respectively .  Maximum  and  minimum  snowfall  was  80.9  in. 
(205.  49  cm)  and  19.7  in.  (  50.04  cm)  in  1953  and  1952 
respectively.  Complete  precipitation  records  for  1971 
through  1974  for  the  Columbia  Icefields  station  located 
approximately  121  km  from  Jasper  show  an  increase  from  78.28 
cm  in  1971  to  115.44  cm  in  1974,  suggesting  that  an  increase 
in  annual  precipitation  has  occurred  since  the  publication 
of  Heusser's  paper. 

During  the  summer  of  1977  when  this  investigation  was 
undertaken  the  weather  was  mainly  cold  and  rainy.  Longer 
periods  of  warm  and  clear  weather  began  during  the  first 
half  of  August  when  most  of  the  data  was  collected.  Figure 
1-2  presents  the  daily  mean  discharge,  average  temperature 
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and  average  per  cent  cloud  for  the  gauging  station 
established  at  the  top  end  of  Beauty  Creek  Flats  for  the 
period  July  21  to  August  15,  1977, 
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1.3.3  Physiography 

The  Sunwapta  River  along  Beauty  Creek  Flats  has  a 
valley  width  of  approximately  0,5  km.  It  is  bordered  on  the 
east  by  Mt.  Sunwapta  and  Tangle  Ridge  at  approximately  3350 
m  and  2900  m  respectively.  On  the  west  are  mountains  of 
slightly  smaller  size.  The  valley  walls  are  heavily  treed 
from  the  level  of  the  gravel  flats  to  the  tree  line  and 
display  avalanche  scars.  The  braided  gravel  flats,  properly 
termed  a  valley  train  (the  equivalent  Icelandic  term  is 
dalsandur  ;  Church,  1972,  p.  4)  are  at  an  approximate 
elevation  of  1550  m.  Mt.  Athabasca  and  the  Columbia 
Icefields  are  visible  to  the  south  and  the  linear  trend  of 
mountains  known  as  the  Endless  Chain  lie  along  the  east  side 
of  the  valley  to  the  north,  out  of  the  study  area. 

1.3.4  Surf icial  And  Bedrock  Geology 

Neither  surficial  nor  bedrock  geology  were  investigated 
in  any  detail  as  they  did  not  relate  to  the  topic  of  the 
investigation.  Inspection  of  aerial  photographs  suggest  that 
the  valley  walls  may  be  veneered  with  lateral  moraine 
remnants  but  this  was  not  corroborated  by  field  work.  Along 
the  south  end  of  Beauty  Creek  Flats  buff  coloured  limestone 
and  shaly  to  silty  limestone  outcrops  were  noted  at  several 
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-F.1fture  1  —  T]^e  daily  mean  discharge,  average  temperature  and  average  percent  cloud  cover  for  the  period 
July  21  to  August  15,  1977  from  the  gauging  station  at  the  top  of  Beauty  Creek  Flats. 


locations  adjacent  to  the  valley  train.  In  one  instance 
silty  limestone  outcropped  as  a  small  island  in  the  gravelly 
alluvium. 

JL  3._5  General  Description 

The  Sunwapta  River  is  proglacial.  Its  meltwaters  are 
derived  principally  from  Athabasca  and  Dome  Glaciers  located 
approximately  20  km  south  of  Beauty  Creek  Flats.  They  are 
outlet  glaciers  of  the  Columbia  Icefields.  Within  the  upper 
20  km  the  Sunwapta  River  receives  additional  discharge 
contributions  from  Stutfield  Glacier  and  an  unnamed  creek 
joining  the  river  from  the  west  and  from  Tangle  Creek  and 
Beauty  Creek  which  join  the  river  from  the  east.  The  former 
two  contribute  meltwater  while  the  water  of  the  latter  two 
is  derived  from  snowmelt  and  pluvial  runoff.  The  unnamed 
creek  and  Beauty  Creek  join  the  Sunwapta  within  the  study 
area  (see  figure  1-3,  p-  12)- 

A  meandering  regime  exists  on  the  braided  channel 
network  of  the  gravel  flats-  As  a  conseguence  of  this,  the 
area  of  active  channels  shifts  from  west  to  east  from  the 
upper  (south)  to  the  middle  and  lower  (north)  reaches 
respectively-  The  areas  of  inactive  channels  adjacent  to  the 
active  channel  areas  have  been  partially  stabilized  by 
grassy  vegetation.  Small  patches  of  stunted  or  infant 
softwood  occur  in  some  of  the  stanilized  areas. 

Three  transects  were  surveyed  on  Beauty  Creek  Flats  in 
order  to  delineate  the  various  morphological  units  of  the 
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valley  train.  The  transects  were  located  at  1.3r  3.3  and  6.4 
km  from  the  southern  end  of  Beauty  Creek  Flats  and  are  shown 
on  figure  1-4  A,  B  and  C  on  p.  14  and  15.  The  morphological 
units  have  been  given  an  identifying  number,  keyed  to  the 
following  scheme: 

1  -  active  channel  area 

2  -  stabilized  area  of  inactive  channels 

3  -  active  channel 

4  -  braid  bar,  major  or  minor 

5  -  non-stabilize d  area  of  discontinuous  channels  and 
standing  water 

6  -  interchannel  area 

The  paucity  of  active  channels  is  a  function  of  the 
time  of  the  year  the  transects  were  surveyed  (early  June) . 

As  observed  in  the  field  the  morphological  units  are 
gradational  with  each  other. 

A  longitudinal  profile  was  established  by  levelling  the 
water  surface  slope  of  the  dominant  active  channel  (see 
figure  1-5,  p.  16).  The  marked  decrease  in  slope  downstream 


is  attributed  to  an  alluvial  fan  located  at  the  northern  end 
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Figure  1-5.  Water  surface  slope  along  the  dominant  active  channel  on  Beauty  Creek  Flats. 
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of  Beauty  Creek  Flats.  The  fan  constricts  the  river  to  a 
single  channel  and  creates  a  local  base  level  for  the 
channels  of  the  study  area.  Within  approximately  2  km 
upriver  of  the  fan  (below  the  lower  reach  of  this  study)  the 
channels  are  bedded  with  coarse  sand  to  fine  gravel  and  from 
eighty  to  ninety  percent  of  the  valley  width  consists  of 
grassy  stabilized  alluvium  (see  figure  1-3).  As  a 
consequence  of  this  decrease  in  slope  there  is  a 
corresponding  decrease  in  the  grain  size  of  the  alluvium 
downstream.  The  valley  train  sediment  is  dominated  by  coarse 
to  fine  gravel  with  sand  occurring  principally  as 
interstitial  grains,  less  frequently  as  small  irregular 
pockets  and  lensoid  bodies  and  rarely  as  thin  discontinuous 
strata  within  the  gravel. 

As  shown  on  figure  1-3,  highway  93  (the  Icefields 
Parkway)  runs  across  the  eastern  margin  of  Beauty  Creek 
Flats.  The  highway  serves  as  a  constraint  on  fluvial 
activity  and  restrains  the  amplitude  of  the  migration  of  the 
areas  of  active  channels.  In  addition,  dykes  composed  of 
very  coarse  rip-rap  have  been  placed  at  regular  intervals 
along  the  road  projecting  obliquely  onto  the  valley  train. 
These  prevent  undercutting  of  the  road  by  interrupting  the 
lateral  migration  of  individual  channels-  Inspection  of 
pre-highway  93  aerial  photographs  show  that  its  construction 
has  caused  the  following: 

1.  shifting  of  the  areas  of  active  channels 

2.  a  reduction  in  the  width  of  the  participating  valley 
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train  by  sealing  off  the  eastern  marginal  area  of  Beauty 
Creek  Flats 

An  increase  in  the  total  discharge  being  conducted  by 
the  braid  channels  occurs  at  the  interface  between  the  upper 
and  middle  reaches  where  the  discharge  of  the  unnamed  creek 
enters  the  Sunwapta  River-  This  is  meltwater  derived  from  a 
small  unnamed  valley  glacier  and  is  thus  opague  with  rock 
flour-  A  second  increase  in  the  total  discharge  on  Beauty 
Creek  Flats  occurs  where  Beauty  Creek  joins  with  the 
Sunwapta.  This  defines  the  start  of  the  lower  reach-  In 
contrast  with  the  unnamed  creek,  the  water  of  Beauty  Creek 
is  clear,  being  derived  principally  from  snowmelt  and 
pluvial  runoff. 

Plates  1  through  12  (p.  19  -  24)  show  the  upper,  middle 
and  lower  reaches  of  Beauty  Creek  Flats  and  previously 
mentioned  salient  features  of  the  study  area- 
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Plate  1:  A  view  looking  west  at  the  lowermost  portion  of  the 
upper  reach  in  April,  1977.  Ice  and  snow  is  still 
present  in  some  of  the  channels . 


Plate  2;  A  view  looking  south  up  the  upper  reach  during  a  high 
flow  period.  Note  the  barely  submerged  bar  in  the 
center  foreground. 
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Plate  3:  A  westward  view  at  the  uppermost  portion  of  the  middle 
reach  in  April,  1977.  Note  the  very  low  flow  in  some 
of  the  channels  and  the  inactive,  stabilized  western 
margin  of  the  valley  train. 


Plate  A:  A  westward  view  at  the  uppermost  portion  of  the  middle 
reach  in  early  July,  1977.  The  increase  in  discharge 
is  indicated  by  comparison  with  plate  3.  Flow  is  from 
left  to  right  (northward). 
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Plate  5:  A  northward  view  down  the  middle  reach  in  July,  1977  . 
Highway  93  is  in  the  foreground. 


Plate  6:  A  northward  view  down  the  lower  reach  in  mid-summer, 
1977.  The  photo  was  taken  from  highway  93. 


22 


Plate  7 :  A  view  of  the  gorge  cut  through  silty  limestone  by  the 

unnamed  creek  on  the  western  side  of  Beauty  Creek  Flats. 
This  creek  marks  the  beginning  of  the  middle  reach. 


Plate  8:  A  view  of  discharge  from  Beauty  Creek  mixing  with  the 
flow  on  Beauty  Creek  Flats.  This  marks  the  beginning 
of  the  lower  reach. 
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Plate  9:  Aeolian  loess  transport  (center  of  photo)  on  the  middle 
reach  of  Beauty  Creek  Flats  on  a  gusty  day.  Note  the 
small  breach  in  the  dyke  in  the  foreground. 
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Plate  10;  Avalanche  scars  on  the  western  valley  wall  of  Beauty 
Creek  Flats . 


Plate  11  :  A  southward  view  from  highway  93  of  the  grassy  area 

(right  center  of  photo)  located  below  the  lower  reach 
of  Beauty  Creek  Flats  (see  fig.  1-2).  The  channel  in 
the  foreground  is  bedded  with  coarse  sand  and  fine 
gravel . 


Plate  12  ;  The  center  portion  of  the  treed  alluvial  fan  at  the 
bottom  of  Beauty  Creek  Flats.  The  grassy  area  in  the 
center  of  the  photo  is  an  abandoned  channel  on  the  fan. 


2-  Developments  In  The  Hydraulic  Geometry  Approach  To 

Channel  Change 


2. 1  Introduction 

Since  the  concept  of  hydraulic  geometry  vas  first 
introduced  by  Leopold  and  Maddock  (1953),  it  has  become  one 
of  the  most  intensively  investigated  areas  in  fluvial 
geomorpholcgy  (Thornes,  1977,  p.  91).  The  literature  dealing 
with  it  is  multi-disciplinary  with  contributions  from 
fluvial  geomorphology,  river  engineering  and  sedimentology. 
Ehodes  (1978,  p.  194)  states  that  at  least  587  sets  of 
hydraulic  geometry  exponents  have  been  computed.  No  attempt 
is  made  here  to  provide  a  comprehensive  review  of  this 
literature,  rather  a  discussion  of  important  developments  in 
the  hydraulic  geometry  approach  to  channel  change  that  have 
resulted  from  the  many  publications  is  presented. 

2^2  Two  Schools  Of  Thought 

One  of  the  most  significant  developments  concerning 
hydraulic  geometry  and  channel  change  involves  the  evolution 
of  two  schools  of  thought  concerning  the  manner  in  which  the 
hydraulic  geometry  relations  of  a  particular  river  or 
physiographic  or  climatic  region  are  dealt  with.  Leopold  and 
Haddock's  paper  in  which  hydraulic  geometry  exponents  were 
averaged  for  seven  river  basins  initiated  a  school  of 
thought  which  searched  for  uniformity  in  channel  response 
from  among  different  physiographic  and  climatic  regions. 


25 


26 


Subsequent  papers  such  as  that  by  Leopold  and  Miller  (1956) 
and  by  Wolman  (1955)  served  to  support  this  approach.  As 
Park  (1977,  p.  134)  states,  a  preoccupation  of  subsequent 
studies  has  been  to  show  how  closely  their  exponents  compare 
to  the  average  values  for  the  mid-western  United  States 
determined  by  Leopold  and  Maddock.  In  1964  this  principle 
gained  further  support  in  Leopold,  Wolman  and  Miller's  text 
entitled  Fluvial  Processes  in  Geomorphology  which  has  become 
a  standard  reference. 

A  second  and  opposing  approach  is  represented  in  two 
recent  publications,  by  Rhodes  (1977)  in  which  332  sets  of 
hydraulic  geometry  equations  were  utilized  and  by  Park 
(1977)  wherein  a  global  variation  in  hydraulic  geometry 
equations  is  examined.  Rather  than  emphasize  only  the 
average  relations  for  a  river  or  region  and  in  so  doing  mask 
differences  which  exist  amongst  its  stations  or  streams, 
respectively,  Rhodes  and  Park  underscore  the  variability 
amongst  the  exponents  from  different  stations  on  the  same 
river,  between  rivers,  or  from  different  physiographic  and 
climatic  regions.  They  stress  the  need  to  avoid  averaging 
and  to  start  concentrating  on  providing  an  explanation  for 
the  local  variations  in  the  values  of  the  hydraulic  geometry 
exponents.  In  conjunction  with  this  they  also  stress  the 
provision  of  detailed  descriptions  of  local  conditions  at 
the  time  of  data  acquisition. 

In  this  investigation  the  latter  school  of  thought  is 
followed.  The  emphasis  is  on  providing  an  explanation  of  the 
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downstream  variation  in  channel  response-  This  approach  is 
supported  by  Knighton  (1975)  who  demonstrated  that  short 
term  temporal  changes  in  channel  morphology  can  occur  at  a 
given  channel  cross  section  during  flow  measurements  which 
constitute  a  major  source  of  scatter  in  a  plot  of  a 
hydraulic  geometry  relation  for  this  section. 


2- 3  General 

Recently  there  has  been  criticism  of  the  conceptual 
background  involved  in  the  conventional  hydraulic  geometry 
approach  to  channel  change-  It  was  mentioned  in  section  1.1 
that  the  conventional  approach  of  hydraulic  geometry  studies 
involved  the  use  of  the  power  function  (log-linear  plot)  to 
express  the  relationship  between  the  dependent  and 
independent  variables.  The  ability  of  such  a  straightforward 
approach  to  interpret  and  forecast  channel  change  has  been 
questioned  by  several  authors  (Thornes,  1970;  Richards, 

1973,  1976).  The  assumptions  involved  in  the  conventional 
approach  were  discussed  by  Knighton  (1977,  p.  101)  and  are 
as  follows: 

1.  Discharge  (Q)  is  the  dominant  independent  variable  to 
which  the  dependent  variables  adjust  and  by  which  other 
independent  variables  (  such  as  grain  size  of  bed 
material  )  are  controlled  to  varying  degrees. 

2.  Power  functions  are  suitable  expressions  of  the 
hydraulic  geometry  relationships. 

A  given  set  of  relations  represents  the  mean  condition 
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of  stream  behaviour  at  a  given  cross  section. 

An  excellent  discussion  of  the  conceptual  problems 
involved  in  the  use  of  this  approach  was  given  by  Thornes 
(1977,  p,  91-100). 

Although  the  conventional  hydraulic  geometry  approach 
to  channel  change  Is  usually  limited  in  its  interpretative 
and  predictive  abilities  by  such  conceptual  difficulties  it 
has  been  and  is  still  accepted  as  a  suitable  means  of 
analysing  the  behavioural  characteristics  of  natural  streams 
(Knighton,  1977,  p.  101).  The  applicability  of  these 
difficulties  to  this  study  is  discussed  in  section  4.2.4. 

Additional  salient  papers  concerning  channel  change  are 
those  of  Schumm  (1960,  1968,  1969)  which  emphasize  the  role 
of  bank  and  bed  sediment  properties  in  determining  how  a 
channel  will  respond  to  an  imposed  discharge.  Schumm  has 
also  recently  (1977)  published  a  text  in  which  channel 
changes  are  dealt  with. 

Langbein  (1964)  presented  his  minimum  variance  theory 
in  which  the  most  probable,  conservative  values  of  the 
hydraulic  geometry  exponents  are  statistically  derived- 


3-  Data  Acquisition 


3,1  Equipment  Used 

Levelling  and  the  determination  of  flow  velocities  with 
current  meters  were  the  major  field  operations  in  this 
investigation.  The  former  provided  the  means  of  obtaining 
accurate  information  on  channel  dimensions  (length,  cross 
section  morphology)  as  well  as  being  used  for  the  long 
profile  and  morphology  transects  of  the  valley  train  (refer 
to  figures  1-5  and  1-4,  respectively).  Information  on  flow 
velocity  was  acquired  by  the  use  of  current  meters  attached 
to  wading  rods.  Flow  velocities  were  obtained  both  as 
primary  data  and  for  combination  with  cross  sectional  area 
for  some  of  the  discharge  values  presented  in  this  report. 
Two  types  of  current  meters  were  used,  A  Price  type  AA 
current  meter  was  used  when  depths  exceeded  0.31  m  (1  ft) 
and  a  pygmy  current  meter  was  used  when  the  depth  of  flow 
was  less  than  0.31  m.  Some  of  the  ancilliary  equipment 
consisted  of  stopwatches,  shovels  and  buckets  (for  channel 
sampling) ,  measuring  tapes  and  a  Stevens  type  I  stage 
recorder-  Plates  13  and  14  on  p.  30  illustrates  one  of  the 
current  meters  and  the  level  used. 
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Plate  13:  Field  assistant  with  a  Price  AA  current  meter  at  the 
gauging  station  at  the  top  of  Beauty  Creek  Flats, 
June,  1977. 


Plate  14 ;  A  photo  of  the  level  used  for  the  field  work  taken  on 
the  upper  reach  of  Beauty  Creek  Flats  in  June,  1977. 
Highway  93  is  in  the  background.  The  flow  is  from 
right  to  left. 
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3.  2  Flow  Conditions  At  The  Time  Of  Data  A  eg  uis  it  ion 

The  data  on  channel  morphology  and  hydraulics  for  the 
three  reaches  of  Beauty  Creek  Flats  were  obtained  during  the 
first  two  weeks  of  August,  1977  (Aug.  4-15),  when  the  summer 
peak  discharge  period  occurred.  Since  the  data  were  obtained 
throughout  the  day  the  diurnal  discharge  variation  was 
significant.  In  this  context,  the  time  of  day  (implying 
rising  or  falling  stage)  was  recorded  for  each  set  of 
channel  data  obtained.  During  data  acquisition  the  stage 
varied  from  slightly  below  bankfull  (usually  pre  10  A.K.)  to 
bankfull  (usually  6-9  P.M.)  .  For  this  reason  the  adopted 
dominant  discharge  of  this  study  is  bankfull  discharge.  The 
discharge  in  the  braid  channels  of  Beauty  Creek  Flats  during 
the  preceeding  part  of  the  summer  had  not  reached  bankfull 
conditions  (except  occasional  brief  periods  of  1-3  days) 
although  considerable  variation  in  stage  had  occurred 
associated  with  varying  weather  conditions.  As  has  been 
observed  for  other  proglacial  rivers,  a  prolonged  period  of 
both  clear  and  hot  weather  (5-6  lays)  are  required  before 
maximum  discharges  occur.  For  the  Sunwapta  River  at  Beauty 
Creek  Plats  there  was  a  1.5  -  2.0  day  lag  between  a  change 
in  weather  and  a  noticeable  braid  channel  response. 

Measurements  of  bed  load  transport  were  not  attempted 
for  the  Sunwapta  braid  channels  ;  however  it  was  noted  that 
significant  gravel  transport  did  occur  throughout  the 
majority  of  the  approximately  two  week  period  during  which 
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the  data  were  obtained-  It  was  not  always  possible  to  see 
such  transport  since  the  rock  flour  in  suspension  obscured 
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visibility.  The  transport  of  sediment  was  noted  by  either 
feeling  the  collisions  with  a  pole  held  in  the  water  against 
the  bottom  or  by  hearing  collisions  during  movement. 
Significant  gravel  transport  was  not  noted  at  any  previous 
time  during  the  summer;  however,  transport  of  sand  sized 
detritus  occurred  at  varying  rates  throughout  the  summer 
depending  on  the  weather  controlled  discharge  conditions. 
During  the  data  acquisition  it  could  be  seen  that  the  most 
active  braiding  was  occurring  in  the  middle  reach-  At  this 
time  it  was  also  noted  that  the  major  source  of  coarse 
detritus  was  upstream  bank  collapse..  As  the  interstitial 
fines  are  washed  out,  the  heavier  gravel  collapses  into  the 
flow-  The  downstream  transport  distance  of  coarse  detritus 
supplied  in  this  manner  has  been  determined  by  other 
investigators  to  be  in  the  order  of  hundreds  of  meters  per 
flood  season  (Butler,  1977,  p.  1072). 

Samples  of  the  suspended  load  being  carried  by  the 
channels  during  the  time  of  data  acquisition  were  not 
obtained  ;  however,  it  was  felt  that  the  maximum  suspended 
load  conditions  encountered  occurred  at  this  time. 

Water  temperature  varied  considerably  between  the  upper 
and  lower  reaches  cf  Beauty  Creek  Flats-  The  average 
temperatures  of  the  active  channel  water  in  the  upper, 
middle  and  lower  reaches  were  6.28,  8-05  and  9.20  C°, 
respectively- 
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Unfortunately  the  total  discharge  carried  through  the 
middle  and  lower  reaches  during  this  time  is  not  known  since 
the  volume  of  water  contributed  by  the  unnamed  creek  and 
Beauty  Creek  (see  figure  1-3,  p.  12)  was  not  monitored. 
However,  the  total  discharge  in  the  upper  reach  during  this 
period  can  be  estimated  from  the  records  of  a  gauging 
station  established  at  the  head  of  Beauty  Creek  Flats  where 
the  flow  was  constricted  to  a  single  channel  by  what 
appeared  to  be  a  bedrock  protrusion  from  the  western  valley 
wall.  Plates  15  and  16,  p.  34  and  35  present  different  views 
of  the  gauging  station  and  the  stage  recorder  apparatus. 

The  discharge  entering  the  upper  reach  during  August 

3 

4-15,  1977  varied  between  13-21  m  /s.  These  values  were 
obtained  from  a  stage-discharge  relation  established  for  the 
station  from  eleven  discharge  measurements  acquired  between 
June  16-26,  used  in  combination  with  the  stage  recorder 
chart  of  the  data  acquisition  interval.  The  stage-discharge 
relation,  referred  to  as  the  rating  equation  is 

Q  =  27.9  S  0,91  r2  =  0.97  Eq-  3.1 

3 

where  Q  is  water  discharge  (m  /s)  and  S  (m)  refers  to  river 
stage.  Figure  3-1  presents  the  graph  of  this  relation-  The 
derived  discharge  range  for  the  upper  reach  should  only  be 
treated  as  a  good  estimate  due  to  the  following  limitations 
arising  out  of  the  derivation  of  the  rating  equation. 

1.  The  eleven  discharge  values  used  in  the  derivation  of 
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Plate  15:  A  downstream  view  of  the  gauging  station  established 
at  the  top  of  Beauty  Creek  Flats.  The  stilling  well 
is  shown  in  the  center  of  the  photograph.  The  channel 
is  approximately  12  m  wide. 
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Plate  16:  A  view  of  the  Stevens  stage  recorder  contained 

within  the  box  attached  to  the  top  of  the  stilling 
well. 
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Figure  3-1.  The  suggested  stage  -  discharge  relation  for  the  gauging  station  established  at  the 

top  end  of  Beauty  Creek  Flats,  Sunwapta  River.  The  data  was  acquired  from  June  18-26, 
1977. 
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the  relation  represent  a  mixture  of  rising  and  falling 
stage  measurements.  Equation  3.1  represents  an  estimated 
steady  state  relationship  for  this  data. 

2-  Although  an  assumption  of  no  standing  pools  of  water 
(zero  stage)  at  zero  discharge  would  seem  to  be 
reasonable  due  to  the  coarseness  of  the  bed  material  no 

3 

measurements  below  9  m  /s  were  obtained,  thereby  forcing 
this  assumption. 

3-  The  maximum  discharge  used  in  deriving  the  relation  was 

3 

approximately  16  m  /s.  To  determine  the  higher 
discharges  prescribed  for  the  upper  reach  a  log-linear 
extrapolation  of  the  rating  curve  is  required.  There  is 
no  evidence  to  suggest  that  the  relation  would  behave 
linearly  at  the  higher  stages  associated  with  the  larger 
discharges  and  in  many  instances  rating  curves  are  found 
to  curve  upward  at  these  stages. 

4.  Practicality  forced  the  establishment  of  the  gauging 
station  just  above  a  bend  in  the  channel  (see  plate  15, 
p.  34).  The  backwater  effects  found  at  such  a  location 
would  slightly  amplify  the  stage  measurements  recorded. 

5.  The  Stevens  stage  recorder  was  not  operating  at  the 
gauging  station  until  July  21.  The  stage  measurements 
made  prior  to  this,  those  used  in  the  derivation  of  the 
rating  equation,  were  done  by  eye  using  a  levelling  rod 
and  slight  inaccuracies  may  have  been  introduced. 

A  section  of  the  stage  recorder  chart  covering  the 

period  August  5-10  is  presented  as  figure  3-2,  p.  39.  Two 
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features  of  the  stage  fluctuations  during  this  period  are  to 
be  noted. 

1.  The  period  August  5  to  mid-day  August  7  clearly 
demonstrates  a  pattern  of  very  rapid  rise  in  stage  from 
late  morning  to  late  afternoon  followed  by  a  slow  fall 
in  stage  from  the  day's  maximum  to  the  late  morning  of 
the  subsequent  day.  This  type  of  pattern  is 
characteristic  of  proglacial  rivers  during  periods  of 
clear  hot  weather  and  demonstrates  how  quickly  the 
source  glaciers  can  respond  to  solar  heating. 

2.  At  approximately  12:30  A.M.  on  August  9  a  sudden 
increase  in  stage  occurred  which  lasted  approximately 
2.5  hrs-,  peaking  at  3:00  A.M. ,  whose  maximum  stage 
exceeded  the  previous  day's  maximum.  Although  not  shown 
in  figure  3-2  a  similar  burst  occurred  between  2  and  3 
A.M.  on  the  morning  of  August  12.  These  bursts  occurred 
subsequent  to  several  days  of  clear  hot  weather  and  are 
interpreted  as  representing  the  passage  of  flood  waves 
originating  in  some  sort  of  plumbing  burst  in  the 
upstream  glaciers  which  occurred  during  the  previous  day 
of  clear,  hot  weather. 

Although  the  unnamed  creek  and  Beauty  Creek  were  not 
monitored  visual  comparison  of  their  flow  magnitudes  to  that 
of  the  gauging  station  during  August  4-15  suggested  that 
their  discharges  individually,  may  approximate  that  of  the 
gauging  station.  This  being  the  case,  a  guess  at  the  range 
of  total  discharge  on  Beauty  Creek  Flats  during  this  period 
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Figure  3-2.  The  stage  recorder  chart  for  the  period  August  5-10, 

1977.  Daily  maximum  and  minimum  stage  is  indicated 
by  dots  on  the  graph.  The  time  scale  increments  are 
6  hrs . 
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would  be  39-63  m  /s,  assuming  that  infiltration  on  the 
valley  train  was  not  significant.  This  assumption  is 
considered  to  be  reasonable.  Calculations  of  Church  (1972, 
P-25)  for  the  outwash  detritus  of  Ekalugad  Fiord  show  the 
groundwater  discharge  at  this  site  to  be  approximately  0.4% 
or  less  of  the  mean  discharge  across  the  gravel  flats. 


3^3  Methodology 

To  facilitate  the  realization  of  the  first  two 
objectives  outlined  on  p.  6,  Beauty  Creek  Flats  was 
separated  into  an  upper,  middle  and  lower  reach  from  each  of 
which  data  on  channel  morphology  and  flow  conditions  would 
be  obtained  and  then  compared.  In  so  doing,  the  first 
objective  is  realized  by  simply  combining  the  morphology  and 
flow  data  of  the  three  reaches.  As  previously  indicated,  the 
upper  reach  is  defined  as  extending  from  the  gauging  station 
location  up  to  but  not  including  the  point  at  which  the 
unnamed  creek  enters  the  valley  train,  a  distance  of 
approximately  2.0  km  (see  figure  1-3,  p.  12).  The  middle 
reach  extended  from  this  position  to  a  position  located  at 
approximately  5.73  km  from  the  beginning  of  Beauty  Creek 
Flats,  a  location  just  above  the  confluence  of  Beauty  Creek 
with  the  valley  train.  The  lower  reach  extends  from  this 
point  to  a  point  at  7. 19  km  down  the  valley  train.  As 
mentioned  previously,  the  last  approximately  2  km  of  the 
valley  train  were  not  included  in  the  lower  reach  since  the 
channels  of  this  area  are  predominantly  sand  bedded  and  have 
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vegetated,  stabilized  banks,  which  would  cause  their 
behavioural  characteristics  to  contrast  sharply  with  those 
of  the  gravelly,  non-cohesive  bed  and  bank  material  of  the 
channels  chosen  from  the  lower,  middle  and  upper  reaches. 

The  three  reaches  of  Beauty  Creek  Flats  are  basically 
similar  with  each  containing  areas  of  active  as  well  as 
inactive  channels  (see  figure  1-4  A,  B  and  C,  p.  1‘  and  15). 
The  principle  difference  in  terms  of  acquiring  channel  data 
is  the  number  of  active  channels  encountered  in  each  reach. 
The  number  of  active  channels  in  each  reach  during  the  data 
acquisition  period  was  never  actually  counted  ;  however  a 
marked  increase  occurred  in  passing  from  the  upper  to  the 
middle  reach  with  only  a  slight  increase  occurring  in  the 
lower  reach  relative  to  the  middle  reach.  Plates  2,  5  and  6 
(p.  19  and  21)  should  be  referred  to  for  a  visual  impression 
of  each  reach. 

In  the  context  of  the  objectives  presented  on  p.  6,  it 
was  desirable  to  obtain  as  much  information  on  channel 
morphology  and  hydraulics  as  manpower  and  time  would  permit. 
In  order  to  accomplish  this  two  groups  of  data  were 
collected,  each  of  which  contained  morphologic  and  hydraulic 
information  on  the  active  channels  of  the  upper,  middle  and 
lower  reach.  The  two  groups  are  primarily  distinguished  by 
the  manner  in  which  their  discharge  values  were  determined. 
Discharges  of  the  first  group  were  determined  using  current 
meter  flow  velocity  measurements  acquired  in  the  field.  This 
group  of  data  is  referred  to  as  the  measured  discharge  data 
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set.  Discharges  for  the  other  group  of  channel  data  were 
determined  theoretically  using  morphologic  and  bed  material 
grain  size  information  from  the  channels.  This  group  is 
referred  to  as  the  theoretical  discharge  data  set.  The 
relation  used  for  the  discharge  calculation  is  presented  in 
the  subsequent  section.  It  is  to  be  noted  that  the  channels 
of  each  group  are  unique;  ie.  the  same  channel  was  not 
investigated  twice.  The  main  reason  for  separating  the  field 
data  into  two  groups  was  the  amount  of  time  consumed  in 
making  numerous  current  meter  discharge  measurements  with 
only  two  people  involved.  For  each  group  the  channels 
selected  from  the  upper,  middle  and  lower  reaches  had  to 
meet  the  following  criteria  to  ensure  that  their  adjustments 
to  the  imposed  discharge  were  representative  of  unrestricted 
channel  change  in  coarse,  non-cohesive  bed  and  bank 
a lluvium. 

1.  A  Well  Defined  Channelized  Flow 

There  were  many  occurrences  of  vaguely  linear  sheets  of 
water  flowing  on  sub- aerial  bar  tops  and  in  interchannel 
areas  which  were  rejected  since  there  was  no  evidence  of 
channelization  associated  with  them. 

2.  Unrestricted  Lateral  Migration  Potential 

Active  channels  which  flowed  for  all  or  part  of  their 
length  against  the  roadside  dykes  or  in  areas  in  which 
the  bank  material  is  to  varying  degrees 


stabilized  by 
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vegetation  were  not  included  in  the  in  vest igation- 
In  order  to  secure  a  representative  sample  a  wide  range 
of  channel  sizes  was  included  in  each  group  of  data.  To 
obtain  the  morphologic  and  hydraulic  information  a  data 
acquisition  format  was  established  for  each  group  of  data 
and  applied  to  the  channels  chosen  for  inclusion  in  one  of 
the  two  groups.  The  format  for  each  group  is  presented 
below. 


Measured  Discharge  Data  Set 

Thirty  channels,  ten  from  each  of  the  upper,  middle  and 
lower  reaches  were  randomly  selected  and  investigated  in  the 
following  manner: 

1.  Once  a  suitable  channel  was  located  a  30.5  m  (100  ft) 
stretch  of  the  channel  was  measured  off  and  marked  by 
placing  wooden  stakes  in  the  gravel  on  the  banks 
adjacent  to  the  start  and  finish. 

2.  Using  driftwood  floats  and  a  stopwatch,  three  thalweg 
surface  velocities  were  determined  by  timing  the  floats 
placed  in  the  thalweg  over  the  30.5  m  distance. 

3.  At  0,  15.2  and  30.5  m  thalweg  velocity  at  0.6  of  the 
flow  depth  (from  the  surface  down)  was  recorded  through 
the  use  of  either  the  AA  or  pygmy  current  meter.  If 
depth  exceeded  0.46  m  (1-5  ft)  the  velocity  at  0.2  and 
0.8  of  the  depth  was  obtained  and  averaged. 

A  velocity-area  (discharge)  measurement  was  obtained 


4. 


44 


from  a  flow  cross  section  located  at  the  midpoint  of  the 
channel  length  (length  was  defined  by  either  a 
bifurcation  and/or  a  confluence).  This  was  achieved  by 
staking  a  rope  across  the  flow,  measuring  the  width  and 
using  loosely  tied  short  lengths  of  flagging  tape  to 
indicate  the  position  of  the  verticals  at  which  current 
meter  velocity  measurements  were  taken.  The  spacing  of 
the  verticals  was  dictated  by  channel  width  and  by 
attempting  to  follow  a  well  established  rule  of  thumb 
which  suggests  that  accuracy  is  maintained  if  the 
discharge  between  adjacent  verticals  is  kept  less  than 
10%  of  the  total  discharge  of  the  cross  section. 

Once  the  velocity-area  measurement  was  completed  the 
equipment  was  gathered  up  and  another  channel  was  sought 
where  the  procedure  would  be  repeated.  For  each  channel  the 
float  velocities  of  step  2  were  averaged  as  were  the  current 
meter  velocities  of  step  3,  giving  an  average  thalweg 
surface  velocity  and  an  average  thalweg  velocity  at  0.6  of 
the  flow  depth,  respectively.  These  two  averages  were  then 
combined  to  give  an  average  velocity  at  0.6  of  the  flow 
depth  to  average  surface  velocity  ratio.  This  ratio  was 
averaged  (0.78)  for  the  thirty  channels  of  this  data  set  and 
used  to  calculate  the  average  velocity  at  0.6  of  the  flow 
depth  for  the  channels  of  the  theoretical  discharge  data  set 
in  which  only  float  velocities  were  determined  (see  below) . 
Also,  the  average  of  the  ratio  of  average  surface  velocity 
to  average  cross  section  velocity  was  determined  (1.67)  for 
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the  calculation  of  the  average  cross  section  velocities  of 
the  theoretical  discharge  data  set  which  are  used  in  the 
determination  of  the  value  of  Mannings  n  (roughness) . 

Theoretical  Discharge  Data  Set 

Thirty  channels,  ten  from  each  of  the  three  reaches 
were  chosen  for  this  group-  Father  than  picking  the  channels 
at  random  from  throughout  a  reach  as  was  done  for  the 
channels  of  the  measured  discharge  group,  seven  transects 
(numbers  1  and  2  in  the  upper  reach,  numbers  3,  4  and  5  in 
the  middle  reach  and  numbers  6  ana  7  in  the  lower  reach  ; 
see  figure  1-3)  were  walked  across  the  valley  train.  Along 
each  of  these,  intersected,  qualifying  active  channels  were 
investigated  in  the  following  fashion: 

1.  Using  the  level  and  levelling  rod  the  channel  length, 
water  surface  slope  and  bed  slope  were  determined  along 
the  channel  centerline. 

2.  An  upper,  middle  and  lower  flow  cross  section  were 
levelled  (the  middle  cross  section  was  not  always 
obtainable) . 

3.  The  water  temperature  was  recorded. 

4.  At  or  near  the  thalweg  of  the  middle  cross  section  a 
grab  sample  of  the  bed  material  was  obtained  with  a 
shovel  and  bucket.  Loss  of  a  large  amount  of  the  fines 
was  unavoidable  with  this  technique.  The  effect  on  the 
grain  size  distribution  is  discussed  in  the  next 
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chapter. 

5-  The  average  thalweg  surface  velocity  was  determined  by 
averaging  the  times  obtained  from  placing  five  floats  in 
the  thalweg  and  allowing  them  to  travel  the  length  of 
the  channel. 

As  mentioned  earlier,  the  discharge  of  each  of  these 
thirty  channels  was  determined  theoretically  using  the 
morphologic  and  grain  size  information  obtained  in  the  above 
manner  in  a  relation  to  be  described  in  detail  in  the 
subsequent  section.  The  average  velocity  at  0.6  of  the  flow 
depth  was  obtained  for  each  channel  by  applying  the  average 
velocity  ratio  determined  for  the  channels  of  the  measured 
discharge  data  set  to  the  average  thalweg  surface 
velocities. 

The  channel  data  obtained  in  the  three  reaches  of 
Beauty  Creek  Flats  according  to  the  two  formats  described 
above  are  presented  in  tabular  form  as  tables  A-1  and  A-2  on 
p.  130  and  132  respectively,  of  appendix  A.  Included  in 
these  tables  are  some  additional  morphologic  and  hydraulic 
parameters  derived  from  the  field  data.  For  the  measured 
discharge  data  set  (table  A-1)  the  derived  parameters  are 
flow  area,  width  to  depth  ratio,  Froude  number  and  channel 
Reynolds  number.  The  derived  parameters  of  the  theoretical 
discharge  data  set  (table  A-2)  are  flow  area,  average  width, 
average  depth,  average  width  to  depth  ratio,  average 
velocity  at  0.6  of  the  flow  depth,  average  bottom  shear 
stress,  Froude  number,  channel  Reynolds  number,  Manning*s 


total  roughness  parameter,  n,  and  of  course  water  discharge. 
The  standard  equations  for  the  derivation  of  the  hydraulic 
parameters  are  presented  in  appendix  A. 

3.4  The  Theoretical  Discharge  Relation 

The  relation  used  to  determine  the  discharge  for  the 
thirty  channels  of  the  theoretical  discharge  data  set  was 
derived  and  its  use  suggested  by  Dr.  G.  Parker  of  the 
Department  of  Civil  Engineering.  The  equation  uses 
dimensionless  depth  and  width  along  with  two  other 
parameters  (see  below)  to  provide  a  dimensionless  discharge. 
From  the  dimensionless  discharge  a  dimensional  discharge  is 
obtained  by  another  relation  (eg.  3.3).  The  equation  is 


~  &  1  7f)  n  cn  * 

Q  =  4.97  d  S B  (1-2^23)  Eq.  3.2 

c  a 

*  d 
7C  Q 

where  d  =  —  =  dimensionless  corrected  depth 

c  D50 
* 

B  =  =  dimensionless  width 

W  =  channel  width 

d  =  corrected  average  depth  for  the  flow  cross  section  ; 

p 

explained  further  below 

S  =  channel  gradient  ;  water  surface  slope  (  S  )  was  used 

w 

a  =  aspect  ratio  =  W/d^ 

Q  =  dimensionless  discharge 
til 

D^.q  =  50  percentile  of  the  grain  size  distribution 
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In  order  to  determine  the  dimensional  discharge  for  the 
hydraulic  geometry  relations  the  following  equation  was  used 
(suggested  by  G.  Parker). 


<  R  8  D50  D50 


Eg.  3.3 


where  Q  =  water  discharge 
R 


(  ps  -  pf  ) 


=  1.70  =  apparent  density  of  the  sediment 


with  p  =  the  sediment  density  (  limestone  density  of  2.7 
gm/cc  was  used  )  and  =  the  fluid  density 


g  =  gravitational  constant 


Equation  3.2  is  based  on  Keulegan's  1938  logarithmic 
velocity  distribution  equation  for  turbulent  open  channel 
flow.  This  eguation  provided  the  average  velocity  for  the 
entire  flow  cross  section  by  approximating  the  cross 
sectional  shape  of  the  channel  as  rectangular.  Parker1 s 
relation  is  essentially  addressing  the  natural  situation  of 
sloping  channel  margins-  Although  the  corrected  depth  (dc ) 
of  the  relation  is  a  theoretical  allowance  for  sloping 
rather  than  vertical  channel  margins,  the  relation  is  still 
not  completely  compatible  with  natural  flow  cross  sections 
since  it  assumes  depth  to  fall  at  the  midpoint  of  the 
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channel  width  (G.  Parker,  personal  communication).  The  flow 

cross  sections  of  the  thirty  theoretical  discharge  channels 

presented  in  appendix  B  clearly  demonstrate  that  the  average 

depth  falls  nearer  to  the  deeper  outside  margin.  The  average 

depth  value  shown  in  table  A-2  of  appendix  A,  represents  a 

channel  average  derived  from  the  average  depths  of  the 

measured  flow  cross  sections  located  at  the  upper,  middle 

and  lower  channel  positions.  Other  than  placing  the  average 

depth  nearer  to  one  margin  rather  than  at  the  midpoint  of 

the  channel  width  these  averages  satisfy  the  definition  of 

d  since  the  distance  down  the  channel  margins  were  measured 
c 

at  each  flew  cross  section  (thereby  incorporating  sloping 
bank  data  into  the  calculations).  For  this  reason  the 
average  depths  of  table  A-2  may  be  considered  as  an 
approximation  of  dc  (suggested  by  G.  Parker).  These  are  then 
adjusted  by  chopping  3.5  of  the  channel  depth  off  of  either 
side  of  the  cross  section  and  recalculating  the  average 
depth  (suggested  by  G.  Parker) .  This  has  the  effect  of 
moving  the  position  of  the  average  depth  for  the  flow  cross 
section  nearer  to  the  midpoint  of  the  channel  width  and 
thereby  satisfying  the  definition  of  corrected  depth  as  best 
as  possible-  The  depths  arrived  at  by  this  procedure  are 
slightly  greater  than  the  field  data  averages  shown  in  table 
A-2.  Since  two  and  sometimes  three  flow  cross  sections  were 
measured  per  channel  for  the  theoretical  discharge  data  set 
the  adjustment  was  made  for  each  flow  cross  section  and  then 
averaged  to  arrive  at  the  estimated  corrected  depth  value 
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for  each  channel. 

To  summarize,  the  determination  of  dimensionless 

discharge  for  the  channels  of  the  theoretical  discharge  data 

set  through  the  use  of  eqs.  3.2  and  3.3  requires  the  use  of 

a  corrected  depth  value  in  equation  3.2.  The  discussion 

above  explains  how  to  adjust  the  average  depth  values  of  the 

field  data  to  arrive  at  an  estimate  of  d  . 

c 

It  has  been  stated  that  equation  3.2  is  based  on 
Keulegan's  logarithmic  velocity  profile  relation-  In  so 
doing  one  is  assuming  that  a  logarithmic  velocity  profile  is 
strictly  applicable  to  the  channels  of  the  theoretical 
discharge  data  set.  An  indication  of  how  applicable  this 
velocity  distribution  is  to  these  channels  may  be  determined 
by  calculating  the  theoretical  surface  velocity  of  the 
thalweg  using  Keulegan's  logarithmic  velocity  distribution 
relation  and  comparing  this  tc  the  empirical  average  surface 
velocity  of  the  thalweg  contained  in  table  A-2  of  appendix  A 
(suggested  by  5.  Parker).  Keulegan's  relation,  when  used  to 
determine  the  surface  velocity  is 

V  =  2.5  u .  In  30  d  Eg-  3-4 

cf  *  — ■— 

D 

90 

where  V  =  theoretical  surface  velocity 

st 

k 

u  =  shear  velocity  =_(  x/p  )  where  x  is  shear  stress  and 
p  is  fluid  density  ;  x  =  average  bottom  shear  stress  is 
substituted  for  x 

d  =  depth  of  flow 

D  Q  =  90  percentile  of  the  grain  size  distribution 
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This  technique  was  tried  on  several  of  the  channels  and 
it  was  found  that  there  was  a  significant  departure  from  a 
logarithmic  velocity  profile.  As  a  consequence,  since 
equation  3.2  was  based  on  Keulegan* s  relation  the  thirty 
theoretically  determined  discharges  required  adjustment  by  a 
correction  factor.  Since  nearly  all  of  the  theoretical 
surface  velocities  were  larger  than  their  empirical 
counterparts  the  correction  factor  defined  as  (suggested  by 
G.  Parker) 


X  =  V  /  V  Eq.  3.  5 

s  st 

where  V  =  measured  surface  velocity  from  table  A-2 

s 

X  =  theoretical  discharge  correction  factor 

was  usually  less  than  unity. 

The  value  of  parameters  used  in  equations  for  the 
calculation  and  correction  of  discharge  values  of  the 
theoretical  discharge  data  set  (refer  to  equations  3.2  to 
3,5)  are  presented  in  table  A-3  of  appendix  A. 

3.5  Laboratory  Procedures 

Laboratory  analysis  was  limited  to  the  sieving  of  the 
grab  samples  of  the  bed  material  acquired  for  the  channels 
of  the  theoretical  discharge  data  set  in  order  to  determine 
the  grain  size  distribution  for  each  sample.  The  following 
procedure  was  used  for  each  of  the  channel  grab  samples. 
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1.  The  entire  sample  (gravel,  sard  and  fines)  is  weighed 
and  then  split  into  a  gravel  and  sand  portion  for 
further  analysis. 

2.  The  gravel  portion  of  the  sample  is  washed  to  remove  any 
adhering  sand,  silt  or  clay  and  is  then  dried  and 
weighed.  The  fines  from  the  wash  were  retained  and 
passed  (wet  sieved)  through  a  230  mesh  sieve  to  extract 
any  sand  sized  material.  The  sand  caught  in  this  manner 
was  added  to  the  rest  of  the  sand  from  the  sample  and 
the  remainder  of  the  wash  material  was  discarded. 

3-  The  sand  fraction  of  the  sample  is  wet  sieved  through 
the  230  mesh  sieve  to  remove  all  silt  and  clay  sized 
material.  The  sand  is  then  dried  and  weighed.  Again  the 
silt  and  clay  in  the  wash  bucket  were  discarded. 

4.  The  percentage  of  the  total  sample  weight  which  was  due 
to  the  silt  and  clay  sized  detritus  is  determined  by 
summing  the  weights  of  the  clean  sand  and  gravel  from 
steps  2  and  3  and  subtracting  from  the  total  sample 
weight  of  step  1 . 

5.  The  gravel  portion  of  the  sample  is  separated  into 
seventeen  size  classes  by  hand  sieving,  the  sediment 
obtained  for  each  class  being  separately  weighed  and 
retained.  The  grain  size  data  is  presented  in  appendix 


C. 

6.  The  sand  portion  of  the  sample  is  broken  down  to  an 
approximately  115  gm  representative  sample,  the  split 
factor  obtained  and  then  separated  into  ten  size  classes 


using  a  Ro-Tap  sediment  shaker.  The  sediment  of  each 
class  is  separately  weighed  and  retained.  This  data  i 
also  presented  in  appendix  C. 

The  grain  size  data  obtained  in  the  above  manner  was 
used  with  an  A.P.L.  interactive  computer  plotting  package 
known  as  8  plot  to  construct  the  cumulative  grain  size 
curves  for  each  of  the  channel  samples.  These  curves  are 
presented  and  discussed  in  the  next  chapter. 


4.  The  Active  Channels  Of  Beauty  Creek  Flats 


4.1  Morphologic  And  Hydraulic  Characteristics  Of  The  Active 
Channels 

The  first  objective  of  the  investigation  as  stated  in 
section  1.2,  p.  6  is  addressed  in  this  section-  Bather  than 
discuss  the  variation  in  each  channel  parameter  over  the 
valley  train  a  statistical  summary  table  (table  4-1, 
p-  56)  is  presented  which  is  derived  from  tables  A-1  and  A-2 
of  appendix  A.  This  table  provides  the  maximum,  minimum, 
mean  and  standard  deviation  of  each  parameter  for  Beauty 
Creek  Flats  as  a  whole,  and  the  mean  and  standard  deviation 
of  each  parameter  for  each  of  the  three  reaches.  Seven 
morphologic  and  seven  hydraulic  parameters  are  described  in 
table  4-1.  The  statistics  on  bed  material  grain  size 
(D  and  D  ,  sediment  ratios,  roughness,  average  bottom 
shear  stress,  slope  and  channel  length  were  determined  using 
only  channel  data  from  the  theoretical  discharge  data  set. 
The  statistics  on  the  remainder  of  the  parameters  contained 
in  the  table  were  determined  by  combining  the  two  data  sets. 
It  should  be  noted  that  in  so  doing  the  statistics  on 
parameters  such  as  width  for  example,  are  determined  by 
combining  channel  average  values  of  the  theoretical 
discharge  data  set  (two  or  three  flow  cross  sections  per 
channel)  with  values  of  the  measured  discharge  data  set 
(only  one  flow  cross  section  per  channel) -  It  is  estimated 
that  approximately  705*  of  the  qualifying  active  channels 


54 


55 


were  investigated  from  each  of  the  three  reaches  and  that 
therefore,  approximately  10%  of  the  qualifying  active 
channels  on  Beauty  Creek  Flats  during  the  period  of  August 
4-15,  1977  are  represented  by  table  4-1, 

•‘•able  4  —  1  shows  that  the  mean  value  of  the  parameters 
progressively  decreases  in  value  from  the  upper  to  the  lower 
reach  of  Beauty  Creek  Flats  with  the  exception  of  Froude 
number,  width  to  depth  ratio,  depth,  width  and  area.  Of 
these,  width  to  depth  ratio  progressively  increases  down 
Beauty  Creek  Flats  while  the  remainder  show  no  consistent 
down  valley  behaviour.  Seven  features  of  the  data  warrant 
discussion. 

1.  Although  the  average  water  surface  slopes  of  the  three 
reaches  are  equal  to  the  corresponding  average  bed 
slopes,  reference  to  table  A-2  demonstrates  that  in  the 
ma jo rit y  of  channels  of  this  data  set  the  former  was 
slightly  greater  than  the  latter. 

2.  A  dominant  control  by  slope  (which  is  itself  determined 
by  the  alluvial  fan  at  the  bottom  of  Beauty  Creek  Flats) 
is  manifest  in  the  fact  that  velocity  decreases  down  the 
study  area  even  though  roughness  (Mannings  n)  is  also 
decreasing.  An  approximate  order  of  magnitude  difference 
exists  between  the  of  the  upper  and  middle  reaches 
and  the  of  the  lower  reach  due  to  the  latter's 
proximity  to  the  alluvial  fan. 

3.  The  mean  values  for  discharge  (Q)  represent  the  average 
discharge  per  channel  in  each  of  the  three  reaches.  As 
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mentioned  above  these  decrease  downstream.  This  occurs 
despite  an  increase  in  the  total  discharge  in  the  middle 
and  lower  reaches  due  to  the  discharge  additions  from 
the  unnamed  creek  and  Beauty  Creek,  respectively. 
Although  the  number  of  active  channels  in  each  reach  was 
never  tallied  this  is  indirect  proof  of  a  progressive 
increase  in  the  number  of  active  channels  down  Beauty 
Creek  Flats. 

A  strong  progressive  decrease  in  grain  size  is  noted  for 

both  the  D 50  and  D  percentiles.  In  both  cases  the 

average  in  the  lower  reach  is  only  approximately  one 

third  of  the  average  in  the  upper  reach.  On  average 

D  =  2.30  D,n.  Grain  size  must  be  considered  a  major 
90  >)U 

factor  in  the  variation  of  channel  behaviour  among  the 
three  reaches  of  Beauty  Creek  Flats. 

The  mean  widths  indicate  that  there  is  a  decrease  in  the 
middle  reach  and  an  increase  (to  a  maximum)  in  the  lower 
reach.  Initially,  the  figure  for  the  lower  reach  may 
seem  a  bit  confusing  since  the  average  discharge  per 
channel  indicates  that  there  is  an  increase  in  the 
number  of  channels  in  this  reach  which  would  cause  one 
to  expect  a  decrease  in  average  width  as  was  the  case  in 
the  middle  reach.  A  consideration  of  the  magnitude  of 
the  increase  in  the  number  of  channels  between  the  upper 
and  middle  and  middle  and  lower  reaches  plus  a 
consideration  of  the  size  of  the  total  discharge 
increase  in  the  middle  (due  to  the  unnamed  creek)  and 
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lower  reacL  (due  to  Beauty  Creek)  explains  this.  Using 
the  average  discharge  per  channel  as  an  index, 
representing  the  increase  in  the  number  of  channels, 
there  is  an  approximately  35%  increase  in  the  number  of 
active  channels  in  the  middle  reach  and  only  an 
approximate  6%  increase  in  the  lower  reach.  In 
conjunction  with  this,  although  the  two  tributaries  were 
not  gauged  it  was  obvious  that  the  discharge 
contribution  from  Beauty  Creek  exceeded  that  from  the 
unnamed  creek.  The  combined  effect  of  a  smaller  increase 
in  the  number  of  channels  and  a  larger  increase  in  the 
total  discharge  resulted  in  the  lower  reach  channels 
being  wider  on  the  average.  Also,  the  reduced  slope  in 
the  lower  reach  would  contribute  to  this  increased 
width. 

In  the  above  context  it  is  reasonable  to  infer  that 
the  lower  reach  channels,  on  average,  should  experience 
overbank  flows  more  frequently  than  the  channels  of  the 
middle  or  upper  reach.  Since  this  would  supply  more  fine 
interstitial  sediment  to  the  inactive  areas  of  the  lover 
reach  it  is  suggested  that  this  may  account  for  an 
upstream  progression  of  the  grassy  area  located  at  the 
bottom  end  of  this  reach  (see  figure  1-3) .  The  upstream 
growth  of  this  area  was  noted  upon  examination  of  air 
photos  taken  in  different  years. 

During  data  acquisition  the  bankfull  discharges  were 
capable  of  moving  only  the  smaller  gravel  sizes  of  the 
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channels-  The  coarser  material  was  moved  for  a  short 

distance  only  when  it  was  dropped  into  the  flow  from 

adjacent  tanks  via  winnowing  of  interstitial  fines. 

It  is  instructive  to  compare  these  observations  of 

sediment  transport  to  the  predictions  provided  by  Baker 

and  Hitter's  (1975)  critical  shear  stress  diagram.  This 

diagram  is  applicable  to  coarse  bedload  material  and  was 

constructed  by  standardizing  the  method  of  critical 

shear  stress  calculation  for  six  sets  of  field 

competence  data.  In  so  doing  they  hoped  to  remove  much 

of  the  scatter  normally  encountered  in  such  plots  where 

data  derived  by  different  techniques  are  grouped 

together.  The  correlation  coefficient  obtained  for  their 

regression  was  0.92  indicating  a  high  degree  of  success 

(see  their  figure  1 ,  p.  976).  Their  scour  diagram  is 

presented  below  as  figure  4-1.  The  average  bottom  shear 

stress  values  of  the  upper,  middle  and  lower  reach  are 

plotted  against  the  average  value  of  each  reach  (see 

table  4-1  for  this  data).  The  plotting  position  of  the 

three  reaches  is  to  the  right  of  the  critical  shear 

stress  line  by  a  significant  distance.  This  indicates 

that  the  D  of  the  channels  in  all  three  reaches  should 
90  - 

have  been  in  transport.  The  critical  shear  stress  for 
the  average  of  the  upper,  middle  and  lower  reach 

channels  as  predicted  by  Baker  and  Hitter's  diagram  is 

2 

approximately  2.20,  1.20  and  0.30  kg/m  ,  respectively. 
This  is  much  less  than  the  calculated  average  bottom 


62 

shear  stress  for  each  reach.  Assuming  the  observations 
on  coarse  gravel  transport  are  accurate  (at  times  these 
were  difficult  due  to  the  amount  of  rock  flour  and 
coarser  fines  being  transported  in  suspension)  the  data 
of  this  study  indicate  that  tne  critical  shear  stress 
values  as  predicted  from  Baker  and  Ritter's  diagram  are 
too  low.  Although  the  imbricated,  underloose  state  of 
the  channel  bed  gravel  may  be  held  partially  accountable 
for  this  discrepancy  the  fact  that  the  data  used  in  the 
construction  of  the  diagram  included  channels  which  also 
would  have  been  underloose  (Fahnestock's  (1963)  channel 
data)  and  the  very  large  differential  between  the 
predicted  critical  shear  stress  values  and  the 
calculated  average  bottom  shear  stress  values  indicates 
that  the  diagram  may  be  inaccurate  for  armoured,  shallow 
and  fast  gravel  bed  channels. 

The  within  reach  variation  in  bed  material  grain  size 
is  represented  in  figures  4-2,  4-3  and  4-4,  p.  65  -  67  which 
present  the  cumulative  curves  of  the  bed  material  for  the 
channels  of  the  theoretical  discharge  data  set  grouped  into 
the  three  reaches  of  Beauty  Creek  Flats.  It  was  mentioned  in 
chapter  three  that  the  method  of  obtaining  the  samples 
(shovel  and  bucket)  forces  a  bias  in  the  resulting  grain 
size  distribution  in  favour  of  the  coarse  portion  of  the 
sample.  The  strong  positive  skewness  indicated  by  the  curves 
in  figures  4-2,  4-3  and  4-4  reflect  this  sampling  bias.  This 
was  unavoidable  since  the  more  sophisticated  eguipment  for 
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Figure  4-1.  Baker  and  Ritter's  (1975,  fig.  1)  plot  of  particle  size 
against  mean  shear  stress  for  coarse  bedload  material 
transported  in  rivers,  canals  and  flood  channels.  In  their 
figure  the  data  points  were  differentiated  into  six  source 
papers.  The  plotting  position  for  the  average  D^q  of  the 
upper,  middle  and  lower  reach  against  the  average  bottom 
shear  stress  of  each  reach  is  indicated  by  u,  m  and  1. 


sampling  bed  material  load  which  may  have  avoided  this 
problem  were  unavailable.  There  are  several  procedures 
available  for  correcting  sample  bias  arising  from  sampling 
technique  (Kellerhals  and  Bray,  1971  and  Potter,  1979), 
however  these  procedures  are  appropriate  only  when  the 
surface  of  the  channel  bed  is  similar  to  the  underlying 
sediment,  or  in  other  words  when  channel  armouring  does  not 
exist  (Potter,  1979,  p.  22).  However,  since  the  samples  are 
dominantly  gravel  (see  table  4-1)  the  weight  of  the  fines 
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lost  by  the  sampling  method  would  have  only  been  a  small 

percentage  of  the  total  weight.  For  this  reason  the  D  and 

50 

Devalues  obtained  from  the  cumulative  curves  are  considered 
to  be  in  excess  of  the  true  values  by  only  a  small  amount. 
The  D  and  D  values  for  the  channels  of  the  theoretical 
discharge  data  set  are  included  in  table  A-2  of  appendix  A. 
It  can  be  seen  from  figures  4-2  to  4-4  that  considerable 
overlap  exists  for  the  curves  of  the  upper  and  middle  reach 
whereas  the  curves  of  the  lower  reach  are  easily 
distinguished  from  these  curves  by  their  finer  grain  sizes. 

Appendix  B  contains  the  flow  cross  sections  obtained 
for  the  channels  of  the  theoretical  discharge  data  set.  A 
within  reach  and  in  some  instances  a  within  channel 
significant  variation  in  flow  cross  section  shape  is 
indicated  in  these  diagrams.  The  large  width  to  depth  ratios 
characteristic  of  channels  in  coarse,  noncohesive  alluvium 
is  apparent  in  these  cross  sections.  In  many  instances  a 
composite  nature  is  revealed  such  as  in  the  case  of  the  top 
and  bottom  flow  cross  sections  of  channel  1  of  transect  4 
(4/1/T  and  4/1/B) .  For  such  channels  it  is  apparent  from  the 
drastic  decrease  in  depth  that  a  smaller  parent  channel  has 
been  enlarged  by  the  scouring  associated  with  overtopping  of 


the  banks. 
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Figure  4-2.  Cumulative  curves  of  the  grab  samples  from  the 

upper  reach  channels  of  the  theoretical  discharge 
data  set. 


grain  size  (phi  notation) 
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Figure  4-3  .  Cumulative  curves  of  the  grab  samples  from  the 

middle  reach  channels  of  the  theoretical  discharge 
data  set. 
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Figure  4-4 .  Cumulative  curves  of  the  grab  samples  from  the 

lower  reach  channels  of  the  theoretical  discharge 
data  set. 
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4  ._2  The  Hydraulic  Geometry  Of  The  Active  Channels 

4.  2.  1  General 

The  second  objective  of  this  investigation,  as  stated 
in  section  1.2,  was  the  interpretation  of  a  variation  in 
channel  behaviour  through  the  development  of  hydraulic 
geometry  equations.  That  objective  is  discussed  in  this 
section  (4.2).  Power  functions  are  presented  (in  4.2.2) 
which  indicate  that  there  is  a  variation  in  channel 
behaviour  among  the  three  reaches  of  Beauty  Creek  Flats-  An 
interpretation  of  the  reason  for  this  variation  is  suggested 
(in  4.2.3)  with  the  emphasis  on  the  explanation  of  the 
variation  in  the  exponents  of  the  hydraulic  geometry 
equations.  The  need  for  such  an  approach  was  emphasized  by 
Park  (1977,  p.  144)  when  he  stated  that  "...in  most 
published  studies  of  hydraulic  geometry  there  is  little 
attempt  to  account  for  variations  in  hydraulic  geometry 
relationships..."  and  "it  is  clear  that  perhaps  the  greatest 
problem  remains  that  of  accounting  for  observed  variations 
in  the  hydraulic  geometry  relationships".  An  explanation  of 
the  variation  in  the  coefficients  of  the  hydraulic  geometry 
relations  is  not  attempted.  Coefficients  reflect  local 
variations  in  sediment  and  flow  characteristics.  Although 
their  interpretation  would  help  explain  the  variation  in 
channel  behaviour  downstream,  the  detailed  information  that 
is  required  for  their  interpretation  was  not  obtained. 


■ 
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The  hydraulic  geometry  relationships  presented  in 
section  4.2.2  were  derived  by  collating  tables  A-1  and  A-2, 
and  then  breaking  the  composite  data  set  down  into  three 
groups  of  data,  one  for  the  upper,  middle  and  lower  reach  of 
Beauty  Creek  Flats.  In  this  manner,  the  variation  in  channel 
behaviour  downstream  could  be  examined.  As  explained  in 
section  3.3  the  channels  were  chosen  at  random  (measured 
discharge  data  set)  and  by  running  transects  across  the 
valley  train  at  various  locations  along  Beauty  Creek  Flats 
(theoretical  discharge  data  set) .  Channels  in  such  a  fluvial 
environment  (braided  gravelly  outwash)  are  characterized  by 
the  regular  occurrence  of  riffles  and  pools  along  the  flow 
path.  The  former  are  shallower,  higher  velocity  areas  found 
at  the  beginning  of  two  channels  at  a  bifurcation  point.  In 
contrast,  pools  are  areas  of  increased  depth  and  decreased 
velocity  found  at  the  confluence  of  two  channels.  As  would 
be  expected  these  two  areas  differ  significantly  in  their 
hydraulic  characteristics.  Bichards  (1976)  has  found  that 
the  rate  of  increase  of  depth  and  velocity  for  a  riffle  area 
is  greater  and  less  than,  respectively,  corresponding  rates 
of  increase  for  a  pool.  Church  (1572,  table  31,  p.  191) 
determined  the  hydraulic  geometry  exponents  of  width,  depth, 
velocity  and  area  for  riffles  and  pools  of  Baffin  Island 
sandar.  Also,  Knight cn  (1975)  has  found  that  riffles  and 
pools  may  be  distinguished  by  differences  in  the  response  of 
flow  resistance,  slope  and  width  to  an  increase  in 
discharge.  For  these  reasons,  when  a  channel  was  intersected 


s 
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on  a  transect  the  flow  cross  sections  were  not  located  in 
either  of  these  two  areas.  Although  it  was  realized  that  the 
differences  between  riffles  and  pools  are  minimized  at 
bankfull  discharge  (Leopold  and  Langbein,  1962),  it  was  felt 
that  such  remaining  variation  might  at  least  partially  mask 
a  downstream  variation  in  channel  response. 

Since  the  emphasis  of  this  investigation  is  the 
explanation  of  the  downstream  hydraulic  geometry  the 
principle  of  constant  frequency  of  discharge  becomes  of 
paramount  importance.  This  principle  was  explained  in 
section  1.  1.  Ideally,  in  a  single  channel  river  situation 
this  principle  would  be  adhered  to  by  examining  the 
discharge  records  of  each  station  along  the  river  and  simply 
picking  those  discharges  (and  accompanying  cross  sectional 
data)  from  each  station's  records  that  have  the  same 
frequency  of  occurrence.  Any  variation  in  the  hydraulic 
geometry  relations  between  stations  would  then  be 
representing  real  differences  and  not  just  differences 
arising  from  comparing  channel  response  at  different 
discharges.  For  a  braided  river,  adherence  to  this  principle 
is  much  more  difficult  since  it  is  not  practical  to  place  a 
gauging  station  on  each  anabranch.  On  Beauty  Creek  Flats,  by 
grouping  the  data  into  three  reaches  and  developing 
hydraulic  geometry  relations  for  each  reach,  the  reach 
becomes  in  effect  a  station.  It  is  important  therefore,  to 
be  able  to  demonstrate  that  for  each  reach  the  average 
discharge  of  the  period  of  data  acquisition  was 
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approximately  constant.  Also,  since  it  was  impossible  to 
refer  to  gauging  station  records,  the  concern  had  to  be  for 
a  similarity  in  average  d ischarge  magnitude  rather  than  a 
similarity  in  frequency  of  occurrence-  As  has  been  stated, 
the  adopted  dominant  discharge  during  the  period  of  data 
acquisition  for  all  three  reaches  was  from  slightly  below 
to,  in  some  instances,  slightly  above  bankfull  discharge. 
Discharge  varied  due  to  both  diurnal  fluctuations  in  stage 
and  weather  induced  differences  in  maximun  discharge  from 
one  day  to  another.  An  indication  of  the  degree  of 
similarity  in  average  discharge  for  the  period  of  data 
acquisition  for  each  reach  may  be  gained  by  reference  to  the 
stage  recorder  chart  from  the  gauging  station  at  the 
beginning  of  Beauty  Creek  Flats  (see  figure  1-3) -  To  do  so, 
one  must  assume  that  a  variation  in  stage  at  this  station 
would  be  matched  by  similar  variations  in  stage  on  the  two 
tributaries.  This  is  felt  to  be  a  reasonable  assumption 
since  all  three  locations  have  a  weather  controlled 
discharge.  In  this  context  then,  the  variation  in  stage  at 
the  gauging  station  may  be  considered  to  be  representative 
of  the  variations  in  stage  on  Beauty  Creek  Flats  as  a  whole. 
The  technique  used  to  determine  the  average  stage  (which  is 
representative  of  discharge)  for  each  reach's  data 
acquisition  period  was  to  firstly  determine  on  which  days  of 
the  period  August  4-15,  the  data  on  the  channels  of  each 
reach  was  obtained.  Once  this  was  done  the  average  stage  for 
each  day  of  the  August  4-15  interval  was  determined  by 
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averaging  the  minimum  and  maximum  stage  of  each  day  as 
determined  from  the  stage  recorder  chart  of  the  gauging 
station-  To  obtain  the  average  stage  for  a  reach's  period  of 
data  acquisition  it  was  necessary  to  weight  the  daily 
average  stages  according  to  the  number  of  channels 
investigated  on  that  particular  day.  In  this  fashion  a 
weighted  average  stage  was  obtained  for  each  day  of  a  given 
reach's  period  of  data  acquisition-  The  weighted  average 
stage  for  the  reach  was  then  obtained  simply  by  averaging 
the  daily  weighted  average  stages  of  that  reach's  data 
gathering  period-  Table  4-2  gives  each  step  of  the  procedure 
for  each  reach.  The  resulting  weighted  reach  averages  are 
0.16,  0-11  and  0-15  m  for  the  upper,  middle  and  lower 
reaches  respectively.  The  similarity  of  these  figures  is 
considered  to  be  sufficient  to  indicate  that  the  variation 
in  channel  behaviour  between  the  three  reaches  of  Beauty 
Creek  Flats  as  indicated  by  the  hydraulic  geometry  exponents 
to  be  presented  subsequently,  represents  actual  reach 
differences  rather  than  differences  attributable  to  the  time 
at  which  the  data  for  each  reach  was  gathered. 


UPPER  REACH 


73 


wt.  reach  average  (m) 

91*0 

wt.  reach  average  (m) 

0.11 

wt .  reach  average  (m) 

0.15 

/-N 

/““N 

E 

E 

E 

'•w' 

CU 

<U 

CU 

00 

to 

oo 

cO 

sr 

CN 

■<r 

CN 

cO 

in 

sO 

r>. 

sf 

co 

CO 

CN 

so 

4-1 

rH 

rH 

rH 

CN 

4-1 

o 

rH 

o 

O 

CN 

4-1 

rH 

rH 

rH 

rH 

c o 

• 

• 

• 

• 

CO 

• 

• 

• 

• 

• 

CO 

• 

• 

• 

• 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 

o 

o 

o 

• 

• 

• 

> 

> 

> 

cO 

co 

co 

• 

• 

• 

+J 

4-1 

4-1 

3 

3 

M 

G 

u 

O 

o 

o 

4-1 

33 

4-) 

4-1 

o 

u 

CJ 

33 

a 

CO 

m 

m 

o 

o 

<0 

CO 

o 

o 

o 

o 

O 

o 

cO 

o 

o 

o 

o 

14-1 

cn 

CN 

CN 

CO 

s 

4-1 

rH 

CO 

rH 

rH 

<5 

4-1 

CN 

CO 

CN 

CO 

• 

• 

• 

• 

pe 

• 

• 

• 

• 

• 

prN 

• 

• 

• 

• 

4J 

o 

o 

O 

o 

4-1 

o 

o 

o 

o 

o 

2 

4-1 

O 

o 

o 

o 

rC 

w 

30 

43 

GO 

x 

GO 

os 

00 

•H 

Q 

•H 

•H 

<u 

Q 

CU 

§ 

CU 

3 

H 

3 

o 

3 

S 

31 

(U 

cu 

cu 

4-J 

4-J 

4J 

to 

co 

CO 

X 

XI 

X) 

co 

CO 

CO 

rH 

m 

n 

-CT 

sO 

rH 

CN 

sO 

CN 

CN 

00 

rH 

■<T 

SO 

>3" 

so 

<u 

(U 

cu 

c 

c 

(3 

c 

G 

C 

cO 

cO 

cO 

43 

43 

4= 

<j 

CJ 

CJ 

E 

E 

E 

'w' 

(U 

<U 

<U 

00 

nt 

00 

CN 

CN 

00 

in 

r-» 

in 

CN 

o 

oo 

so 

O 

CN 

CO 

m 

r^» 

cO 

'O' 

in 

r>* 

sO 

co 

sO 

in 

SO 

n 

4-1 

• 

• 

• 

• 

4-1 

• 

• 

• 

• 

• 

4-> 

• 

• 

• 

• 

CO 

o 

o 

o 

o 

CO 

o 

o 

o 

o 

o 

CO 

o 

o 

o 

o 

• 

• 

• 

> 

> 

> 

CO 

c0 

CO 

m 

CN 

co 

vO 

r^. 

00 

CO 

Os 

o 

■vT 

in 

<u 

fH 

rH 

cu 

rH 

rH 

(U 

rH 

rH 

rH 

4-1 

• 

4-) 

• 

4-1 

• 

cO 

00 

cO 

00 

CO 

oo 

T3 

3 

X 

0 

X 

3 

< 

<3 

< 

Table  4-2.  The  procedure  used  to  estimate  the  degree  of  similarity  in  flow  magnitude  between  the  data 
gathering  periods  of  the  upper,  middle  and  lower  reach  of  Beauty  Creek  Flats. 
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4.2.2  Presentation  Of  The  H ydraulxc  Geometry  Relations 

In  this  section  the  power  functions  developed  for  the 
upper,  middle  and  lower  reaches  of  Beauty  Creek  Flats  are 
presented.  Four  channel  characteristics  are  regressed 
against  discharge  for  each  reach.  The  comparison  of  the 
exponent  of  the  power  function  for  a  given  parameter,  for 
example  width,  to  corresponding  exponents  in  the  other  two 
reaches  is  an  indication  of  the  magnitude  of  the  variation 
in  channel  response  between  the  three  reaches  with  respect 
to  this  particular  channel  characteristic.  By  combining  this 
with  similar  comparisons  done  for  the  other  channel 
characteristics  the  amount  of  variation  in  channel  behaviour 
down  Beauty  Creek  Flats  is  indicated. 

It  is  important  to  distinguish  between  the  average 
value  of  a  particular  parameter  for  a  reach  as  is  given  in 
table  4-1  and  the  rate  of  change  of  that  parameter  with  an 
increase  in  discharge  which  is  indicated  by  the  exponent  of 
the  power  function.  A  given  parameter  may  increase  or 
decrease  with  an  increasing  discharge  (positive  or  negative 
correlation  respectively) .  A  set  of  channels  or  a  reach  (in 
the  case  of  the  average  values  given  in  table  4-1)  may  be 
shallower  than  another  but  may  have  a  greater  rate  of  depth 
increase  with  increasing  discharge.  In  this  case,  the 
shallower  reach  may  have  the  smallest  average  depth  but  the 
largest  positive  exponent  for  its  power  function- 

Ideally,  when  hydraulic  geometry  relations  are 
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developed  for  a  station  the  stream  data  from  which  they  are 
derived  has  been  obtained  by  repetitive  flow  measurements 
over  a  period  of  time  at  exactly  the  same  flow  cross 
section.  For  the  three  reaches  of  Beauty  Creek  Flats  only 
one  set  of  measurements  was  obtained  per  channel.  Twenty 
channels  were  investigated  in  each  reach  and  combined  to 
constitute  the  data  set  from  which  the  hydraulic  geometry 
relations  for  that  reach  were  derived.  This  approach  in 
effect,  constitutes  taking  stream  measurements  at  various 
flow  cross  sections  in  tne  \icinity  of  a  designated  station. 
This  is  not  an  invalid  approach  but  the  amount  of  scatter 
would  be  expected  to  be  significantly  greater  than  if  the 
data  came  from  the  same  flow  cross  section  (Wolman,  1955,  in 
Fahnestock,  1963,  p.  13). 

As  mentioned  above,  power  functions  were  determined  for 
four  channel  characteristics  for  each  reach.  The  four 
characteristics  were  width,  average  depth,  average  surface 
velocity  and  water  surface  slope.  In  all  cases  the 
independent  variable  was  discharge.  Power  functions  for  the 
same  four  parameters  were  also  developed  for  Beauty  Creek 
Flats  as  a  whole  by  combining  the  data  sets  of  the  three 
reaches. 

A  standard  notation  exists  for  the  coefficients  and 
exponents  of  the  power  functions  lor  the  four  channel 
characteristics  to  be  presented.  The  general  form  of  the 
relation  for  each  of  these  parameters  is  as  follows: 
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b 

W  =  aQ 

f 

d  =  cQ 

m 

V  =  kQ 

z 

S  =  tQ 

In  the  subsequent  interpretation  (section  4.2,3)  the 
exponents  for  a  particular  relationship  may  be  referred  to 
by  their  corresponding  letter. 

The  power  function  relationships  were  derived  from  base 
10  log-log  plots  of  the  data  on  each  reach.  The  use  of 
logarithms  is  convenient  primarily  due  to  the  fact  that 
river  discharge  tends  to  be  log-normally  distributed  (Blatt, 
Middleton  and  Murray,  1972,  table  3-2,  p.  39).  Histograms  of 
both  transformed  and  untransformed  data  were  examined.  The 
logarithmic  transformation  produced  a  significantly  more 
normal  distribution  for  the  discharge  and  width  data  while 
only  minor  improvement  was  noted  in  the  distributions  of  the 
other  parameters.  The  normalizing  effect  of  the 
transformation  makes  the  data  more  amenable  to  various 
statistical  tests  which  require  an  assumption  of  normality- 
The  power  functions  determined  for  each  reach  and  for 
Beauty  Creek  Flats  as  a  whole  are  presented  below  followed 
by  the  corresponding  graphs  as  figure  4-5  A,  B,  C  and  D.  It 
will  be  noted  that  although  twenty  chancels  were 
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investigated  in  each  reach  the  water  surface  slope  relations 
were  determined  with  only  ten  data  points  (only  nine  for  the 
slope  relation  in  the  lover  reach)  .  This  results  from 
acquiring  slope  measurements  for  the  theoretical  discharge 
data  set  only  (see  section  3.3).  It  should  also  be  noted 
that  the  average  depth-discharge  relation  was  determined  by 
regressing  the  average  depth  values  of  tables  A-1  and  A-2 
against  discharge  rather  than  the  corrected  depth  values 
which  were  used  in  the  theoretical  determination  of 
disch  arge. 
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Figure  4-5  A.  Log-log  plots  (4)  of  the  hydraulic  geometry  relations 

developed  for  the  upper  reach  of  Beauty  Creek  Flats; 
the  previous  page  is  included  in  this  figure. 
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Figure  4-5  B.  Log-log  plots  (4)  of  the  hydraulic  geometry  relations 

developed  for  the  middle  reach  of  Beauty  Creek  Flats; 
the  previous  page  is  included  in  this  figure. 
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Figure  4-5  C.  Log-log  plots  (4)  of  the  hydraulic  geometry  relations 

developed  for  the  lower  reach  of  Beauty  Creek  Flats; 
the  previous  page  is  included  in  this  figure. 
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Figure  4-5  D.  Log-log  plots  (A)  of  the  hydraulic  geometry  relations 

developed  for  Beauty  Creek  Flats  as  a  whole;  the 
subsequent  page  is  included  in  this  figure. 
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— —  Interpretation  Of  Downstream  Variation  In  Channel 
Behaviour 

4 . 2 . 3 . 1  Met hodolo  qy 

The  interaction  between  the  morphologic  and  hydraulic 
parameters  that  determine  a  channel's  response  to  a 
variation  in  the  sediment-water  mixture  is  complex  and  as 
Eichards  (1973,  p.  894)  states  "achieving  a  point  of  entry 
into  the  system  for  the  purpose  of  analysis  is  extremely 
difficult".  One  cannot  pick  a  particular  parameter  or 
several  parameters  and  state  that  only  these  will  affect  the 
response  of  a  certain  channel  characteristic  to  a  variation 
in  discharge.  The  parameters  are  thoroughly  interrelated. 

One  must  choose  dominant  or  control  parameters  from  the 
available  data;  those  which  most  adequately  explain  the 
variation  in  that  particular  channel  characteristic. 
Individually,  the  parameters  may  display  opposing  effects  on 
the  channel  characteristic  being  interpreted,  that  is,  some 
may  tend  to  make  it  increase  with  increasing  discharge  while 
others  would  be  expected  to  produce  a  decrease.  It  is  the 
dominant  parameters  which  exert  the  strongest  control  and 
which  therefore  determine  the  net  behaviour  of  that 
particular  channel  cha racter istic. 

In  the  subsequent  section  an  interpretation  is 
suggested  for  the  downstream  variation  in  each  one  of  the 
four  channel  characteristics  that  were  regressed  against 
discharge.  As  such  the  interpretation  in  section  4. 2-3.2  is 
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a  downstream  interpretation  as  opposed  to  an  at-a-station  or 
at~a~station  downstream  interpretation  (the  latter  is 
presented  in  section  4-2- 3- 4).  The  distinction  between  these 
types  of  interpretation  was  provided  in  section  1.1. 

4. 2. 3. 2  Downstream  Hydraulic  Geometry 


Kidth 

re  ach 

upper 

riddle 

lower 


b  value 
0.39 
0.47 
0.  45 


dominant  par ame  te  rs:  ,  D 

Cn  Eeauty  Creek  Flats  the  rate  of  increase  of  width 

with  increasing  discharge  increases  in  the  middle  reach 

channels  and  decreases  slightly  in  the  channels  of  the  lower 

reach.  Despite  the  increase  in  total  discharge  in  the  middle 

reach  due  to  the  discharge  contribution  from  the  unnamed 

creek  the  concurrent  increase  in  the  number  of  active 

channels  results  in  an  average  discharge  per  channel  which 

3  3 

is  less  than  that  in  the  upper  reach  (2.14  m  /s  vs  3.27  m 
/s) .  In  light  of  the  larger  b  value  for  the  middle  reach 
this  would  initially  seem  anomalous,  however  it  can  be 
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explained  by  reference  to  the  t  and  values  for  the 

o  50 

upper  and  middle  reach  channels.  Ihe  fact  that  the  middle 

reach  channels  are  increasing  their  width  at  a  greater  rate 

indicates  that  the  differential  between  the  x,  (bank  shear 

b 

stress)  and  the  x^  (critical  shear  stress)  is  larger  than  it 

is  in  the  upper  reach  despite  the  lower  average  discharge 

-22 
per  channel  and  the  lower  x  (18.67  kg/c  vs  31.69  kg/m  ; 

o 

average  bottom  shear  stress  is  taken  as  representative  of 

average  bank  shear  stress).  This  is  possible  due  to  a  marked 

decrease  in  the  D  (decreased  by  approximately  30%)  in  the 

middle  reach  channels  which  causes  the  banks  to  be  more 

easily  eroded  (it  is  assumed  that  the  downstream  change  in 

the  size  distribution  of  the  channel  bed  is  representative 

of  the  downstream  change  in  the  size  distribution  of  the 

banks) .  In  short,  the  D  decrease  overrides  the  decrease  in 

the  average  discharge  per  channel  and  average  bottom  shear 

stress  (x  )  in  the  middle  reach  producing  a  x,  /x  ratio 
o  be 

which  is  larger  than  it  is  in  the  upper  reach  resulting  in  a 
more  rapid  increase  in  width.  Note  that  this  is  the  case 
even  though  the  average  width  in  the  middle  reach  is  less 
than  that  in  the  upper  reach;  ie.  more  channels  of  smaller 
widths  as  opposed  to  a  few  large  channels. 

In  passing  from  the  middle  to  the  lower  reach  the 
behaviour  of  the  median  grain  size  (  d50)  ,  the  average  bottom 
shear  stress  (x  )  and  the  average  discharge  per  channel  is 
the  same  as  that  which  occurred  in  passing  from  the  upper  to 
the  middle  reach  in  as  far  as  the  sense  of  their  response  is 
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concerned-  However,  the  decrease  in  is  much  greater  (57% 

vs  30%) ,  the  decrease  in  to  is  less  (38%  vs  41%)  and  the 
decrease  in  the  average  discharge  per  channel  is  less  (6%  vs 
35%) .  Rather  than  suggest  a  slight  decrease  in  the  rate  of 
increase  of  channel  width  the  data  of  the  lower  reach  would 
lead  one  to  expect  an  increase  in  the  value  of  b.  This  would 
be  expected  since  neither  the  average  bottom  shear  stress 
(considered  to  be  indicative  of  the  relative  magnitudes  of 
the  bank  shear  in  the  three  reaches)  nor  the  average 
discharge  per  channel  experiences  as  large  a  decrease  as 
occurred  in  the  upper  to  middle  reach  transition  yet  the 
drops  drastically;  ie.  the  erosive  energy  does  not 
drop  as  much  and  the  bank  material  grain  size  drops  by  a 
much  greater  amount  than  that  which  occurred  from  the  upper 
to  middle  r<  .ch,  which  should  make  the  banks  easier  rather 
than  harder  to  erode.  It  is  suggested  that  the  controlling 
factor  is  the  percentage  of  interstitial  sand  in  the  bank 
gravel.  Assuming  the  percentage  of  sand  in  the  channel  grab 
samples  of  the  theoretical  discharge  data  set  is 
representative  of  the  amount  of  sand  in  the  channel  banks, 
the  increase  in  the  percentage  of  sand  from  the  upper  to  the 
middle  reach  is  approximately  30%  (5.41%  vs  7.69%)  and  the 
increase  from  the  middle  to  the  lower  reach  is  approximately 
53%  (7.69%  vs  16.23%).  The  ultimate  reason  for  this  increase 
is  of  course  the  increasing  proximity  to  the  alluvial  fan  at 
the  bottom  of  Beauty  Creek  Flats  which  lowers  the  competence 
of  the  channels  by  decreasing  the  slope.  The  bank  gravel  of 
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the  middle  channels  having  less  interstitial  sand  would  be 
of  a  more  open  work  nature  and  therefore  more  susceptible  to 
gravity  collapse  initiated  by  bank  shear  and  undercutting 
than  the  channels  of  the  lower  reach.  It  is  felt  that  the 
interstitial  sand  serves  a  supportive  function  which  becomes 
a  significant  factor  affecting  a  slight  decrease  in  the  rate 
of  increase  of  width  in  the  lower  reach  simply  by  virtue  of 
its  greater  abundance.  Fahnestock  (1963,  p.  30)  suggests  the 
same  role  for  the  interstitial  sediment.  The  sand  does  not 
add  any  cohesiveness  but  rather  holds  the  gravel  size 
sediment  more  securely  in  place  by  tighter  packing. 

Briefly  stated  then,  the  increased  sand  percentage  more 
than  counterbalances  the  variation  in  the  other  three 
mor phologic/hydraulic  factors  which  occurs  in  the  lower 
reach  and  which  otherwise  would  be  expected  to  increase  the 
rate  of  increase  or  width. 


Average  Depth 


re  ach 

f  value 

upper 

0.28 

middle 

0.25 

lower 

0.36 

d  ominant  parameters : 
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The  average  depth-discharge  relation  for  the  three 
reaches  of  Beauty  Creek  Flats  indicates  that  the  rate  of 
increase  of  average  depth  decreases  slightly  in  the  middle 
reach  channels  and  increases  markedly  in  the  channels  of  the 
lower  reach. 

The  dominant  parameters  controlling  the  variation  in  f 
down  Beauty  Creek  Flats  are  felt  to  be  tq  and  D^q  .  in 
passing  from  the  upper  to  the  middle  reach  the  drops  by 

30%  (mentioned  above) ,  however  the  of  the  middle  reach 
channels  decreases  by  41%.  This  greater  decrease  in  average 
bottom  shear  stress  is  interpreted  as  overriding  the  drop  in 
median  grain  size,  resulting  in  an  f  value  slightly  lower 
than  that  of  the  upper  reach.  In  the  lower  reach  a  smaller 
drop  in  (38%)  and  a  much  larger  decrease  in  (57%) 

occurs  with  the  net  result  being  more  erodible  channel  beds 
and  an  accompanying  increase  in  the  rate  of  increase  of 
channel  depth.  If  the  interpretation  is  to  be  viable  then 
one  must  assume  that  the  increase  in  the  percentage  of  sand 
will  affect  the  width  response  (discussed  previously)  but 
not  that  of  depth.  Such  may  be  the  case  since  the  intensity 
and  consistency  of  winnowing  activity  would  be  greater  on 
the  channel  beds  than  on  the  banks  resulting  in  a  relative 
absence  of  interstitial  fines  in  the  coarse  surficial  layer 
of  channel  bed  sediment.  This  situation  is  entirely 
reasonable  since  the  occurrence  of  armoured  beds  in  gravel 
streams  is  a  very  common  phenomenon.  The  fact  that  the 
channel  bottom  grab  samples  contained  enough  sand  to  make 
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them  useful  for  making  inferences  as  to  the  percentage  of 
sand  in  the  channel  banks  (see  width  interpretation)  does 
not  indicate  an  absence  of  armouring  in  these  channels.  This 
is  a  function  of  the  shovel  and  bucket  sampling  technique 
which  resulted  in  a  penetration  of  the  channel  bed  and  the 
incorporation  of  more  than  just  the  armoured  surface  layer 
into  the  sample. 


Average  Surface  Velocity 


reach 

upper 

middle 

lower 


m  value 
0.  23 


0.  20 

0.15 


dominant  parameter:  S 

The  interpretation  of  the  variation  in  m  down  Beauty 
Creek  Flats  is  straightforward  and  is  an  excellent  example 
of  the  influence  of  a  more  dominant  factor  in  controlling 
the  behaviour  of  a  particular  channel  character istic.  Table 
4-1  indicates  that  the  average  roughness  value  decreased 
down  Beauty  Creek  Flats,  a  not  unexpected  response  in  light 
of  the  decrease  in  .  Therefore,  if  one  were  to  consider 
only  roughness  the  expected  behaviour  of  m  would  be  to 
increase  down  the  study  area.  However,  rather  than  increase 
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from  the  upper  to  the  lower  reach  m  decreases.  This 

behaviour  is  attributed  to  the  decreasing  value  of  the  water 

surface  slope  as  the  alluvial  fan  is  approached.  The  water 

surface  slope  of  the  middle  reach  is  only  slightly  smaller 

than  that  of  the  upper  reach  (0.010  vs  0.011)  and  this  is 

reflected  in  a  marginally  smaller  m  value  (0.20  vs  0.23).  In 

the  lower  reach,  however,  the  S  is  approximately  an  order 

w 

of  magnitude  less  than  it  is  in  the  middle  reach.  This 
decrease  is  similarly  reflected  in  the  m  value  by  a  drop  to 
0.15.  Therefore,  although  velocity  increases  with  increasing 
discharge  in  all  three  reaches  the  rate  of  increase  is  the 
smallest  in  the  lower  reach  due  to  the  proximity  to  the 
alluvial  fan. 

A  contributing  factor  to  the  decrease  in  m  in  the 

middle  reach  is  the  decrease  in  average  channel  depth  of 

this  reach  (see  table  4-1).  The  latter  decrease  is  as  would 

be  predicted  by  the  Manning  eguation  (see  eg.  A.  4  of 

appendix  A).  The  average  depth  remains  approximately 

constant  in  the  lower  reach  and  therefore  would  not 

compliment  S  t  in  producing  a  decrease  in  V  . 

w  s 


Water  Surface  Slope 

The  slope-discharge  relations  for  Beauty  Creek  Flats 
are  very  poor.  This  is  indicated  by  the  very  low  values  of 
the  coefficient  of  determination  presented  in  figure  4-5  A 
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to  D.  For  this  reason  a  detailed  explanation  of  exponent 
variations  is  not  presented-  The  lack  of  significance  for 
these  regressions  indicates  that  slope  is  only  minimally,  if 
st  all,  dependent  on  discharge  and  that  it  is  therefore, 
effectively,  an  independent  parameter. 

Allowing  the  water  surface  slope  to  be  representative 
of  the  bed  slope  (this  is  justifiable  on  the  basis  of  the 
similarity  of  the  two  slopes  shown  in  table  A-2)  the 
following  explanation  of  the  poor  correlation  between  slope 
and  discharge  is  offered. 

In  most  streams  slope  is  readily  adjusted  by  scour 
and/or  deposition  on  the  channel  bed.  Such  changes  usually 
constitute  minor  modifications  by  the  stream  in  an  attempt 
to  establish  an  equilibrium  with  the  sediment- water 
discharge.  On  Beauty  Creek  Flats  the  slope  of  the  gravel 
anabranchs  is  suggested  to  be  in  a  state  of  nonadjustment 
due  to  the  inability  of  the  majority  of  flows  to  move  the 
sediment.  This  inability  may  be  attributed  to  the  coarse 
nature  of  the  sediment  in  combination  with  the  low  areal 
slopes  of  the  three  reaches  due  to  the  establishment  of  a 
base  level  by  the  alluvial  fan  at  the  bottom  of  Eeauty  Creek 
Flats.  The  degree  of  nonadjustment  of  slope  can  be  expected 
to  vary  from  one  anabranch  to  another.  This  lack  of  slope 
adjustment  by  discharge  is  reflected  in  the  large  amount  of 
scatter  shown  on  the  slope- discharge  plots. 

On  Beauty  Creek  Flats  the  slopes  of  the  channels  are 
not  a  function  of  discharge.  The  downstream  decrease  in 
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slope  is  a  function  of  increasing  proximity  to  the  alluvial 
fan-  The  within  reach  and  in  particular  within  transect 
variation  in  slope  is  suggested  to  be  largely  due  to 
variations  in  the  size  distribution  of  the  channel  bed 
sediment  which  are  in  turn  the  result  of  variations  in  the 
hydraulic  characteristics  of  the  channels. 

4 . 2- 3. 3  Summary  Of  The  Downstream  Hydraulic  Geometry 

The  hydraulic  geometry  relationships  developed  for  the 
three  reaches  of  Beauty  Creek  Flats  indicate  a  significant 
variation  in  channel  response  down  the  study  area.  The  rates 
of  change  of  width,  average  depth  and  average  surface 
velocity  were  compared  using  the  exponents  of  the  hydraulic 
geometry  relationships-  The  independent  parameter  in  all 
cases  was  water  discharge. 

Width 

The  ease  of  bank  collapse  in  the  gravelly  alluvium 
resulted  in  the  primary  adjustment  to  an  increasing 
discharge  being  an  increase  in  channel  width.  This  is 
indicated  by  the  higher  values  of  the  b  exponents  in  all 
three  reaches. 

The  rate  of  increase  of  width  increased  from  the  upper 
(0-39)  to  the  middle  reach  (0.47)  and  decreased  slightly 
from  the  middle  to  the  lower  reach  (0-45).  The  reason  for 
this  variation  is  attributed  to  a  dominant  control  by  a 
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decreasing  D  in  the  former  case  and  by  an  increasing 
percentage  of  interstitial  fines  in  the  latter. 

The  fact  that  the  middle  reach  channels  have  the 
largest  b  value  and  that  coarse  material  supplied  by 
upstream  bank  collapse  is  characteristically  transported 
only  several  tens  of  meters  downstream  in  such  channels 
suggests  that  the  highest  rate  of  sediment  deposition  is 
occurring  in  the  middle  reach.  During  the  period  of  data 
acquisition  it  was  obvious  that  the  most  active  braiding  was 
occurring  in  this  reach.  Since  deposition  of  coarse  bedload 
and  the  resultant  deflection  of  the  flow  in  two  directions 
is  a  cause  of  braiding  the  field  observations  support  the 
conclusion  for  the  highest  sediment  deposition  rate  in  the 
middle  reach.  Also  supporting  this  conclusion  is  the  fact 
that  the  lowest  rate  of  increase  of  depth  (0.25)  occurred  in 
the  middle  reach.  The  implication  here  is  that  the  rate  of 
increase  of  width  may  be  used  as  an  indicator  of  the 
intensity  of  braiding  in  the  gravelly  glaciofluvial 
environment. 

In  the  above  context,  it  is  apparent  that  erodible 
banks  are  the  primary  cause  of  the  braided  pattern  on  Beauty 
Creek  Flats. 

Av  erage  Depth 

In  all  three  reaches  the  secondary  adjustment  to  an 
increase  in  discharge  was  by  an  increase  in  average  channel 
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depth  as  the  f  values  indicate.  The  rate  of  increase  in 

average  depth  decreases  from  the  upper  (0-28)  to  the  middle 

reach  (0.25)  and  increases  considerably  in  the  lower  reach 

(0.36) .  The  controlling  parameters  over  this  variation  are 

suggested  to  be  average  bottom  shear  stress  (t  )  and  average 

o 

grain  size  (  D^)  .  The  net  behaviour  of  average  depth  is 
determined  by  the  relative  magnitude  of  the  decrease  in 
these  two  parameters  which  occurs  in  passing  from  the  upper 
to  the  middle  reach  and  from  the  middle  to  the  lower  reach. 
In  the  former  transition  the  decrease  in  average  bottom 
shear  stress  (41%)  exceeds  the  drop  in  D^q  (30%)  producing 
the  smaller  f  value  in  the  middle  reach.  In  the  middle  to 
lower  reach  transition  the  reverse  occurs  and  the  drop  in 
D50  (57%)  exceeds  the  decrease  in  average  bottom  shear  stress 
(38%)  which  allows  the  rate  of  increase  of  average  depth  to 
increase.  The  percentage  of  interstitial  fines  is  not 
considered  important  in  determining  the  depth  response  due 
to  their  removal  from  the  surficial  layer  of  coarse  channel 
bed  sediment  by  winnowing  which  results  in  an  armoured  bed. 

Average  Surface  Velocity; 

In  all  three  reaches  the  tertiary  adjustment  to  an 
increasing  discharge  is  an  increase  in  the  average  surface 
velocity  of  the  channels.  The  rate  of  increase  of  average 
surface  velocity  decreases  from  the  upper  (0.23)  to  the 
middle  reach  (0-20)  and  from  the  middle  to  the  lower  reach 
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(0.15).  The  dominant  influence  over  this  progressive 
decrease  is  the  decrease  in  average  water  surface  slope  down 
Beauty  Creek  Flats  created  by  the  alluvial  fan  at  their 
downstream  end  which  restricts  the  flow  to  one  channel  and 
creates  a  local  base  level  for  the  upstream  channels.  A 
small  difference  between  the  average  slopes  of  the  upper  and 
middle  reaches  (0.011  and  0.010  respectively)  produced  the 
similarly  small  decrease  in  the  m  value  in  the  middle  reach 
and  an  approximate  order  of  magnitude  difference  between  the 
average  water  surface  slopes  of  the  middle  and  lower  reaches 
caused  the  larger  drop  in  the  m  value  in  the  lower  reach. 

The  smaller  values  of  these  velocity  exponents  relative  to 
the  m  values  obtained  in  such  hydraulic  geometry  studies  as 
that  by  Leopold  and  Maddock,  1953  and  Kolman,  1955  whose 
average  m  values  were  0.34  and  0.55  respectively,  is  a 
further  indication  of  the  relative  instability  of  the  Beauty 
Creek  Flats  channels  (see  Rhodes,  1977,  p.  80). 

Slope 


The  very  low  correlation  for  the  slope-discharge  power 
function  relations  developed  for  Eeauty  Creek  Flats 
indicates  that  slope  may  be  considered  as  an  independent 
parameter. 

Slope  is  suggested  to  be  in  a  state  of  nonadjustment 
due  to  the  inability  of  most  of  the  flows  transmitted  by  the 
channels  to  move  the  coarse  gravelly  sediment. 
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The  downstream  decrease  in  slope  is  attributed  to  an 
increase  in  proximity  to  the  alluvial  fan  at  the  bottom  of 
Beauty  Creek  Flats  and  a  within  reach  and  within  transect 
variation  is  suggested  to  be  mainly  due  to  variations  in  the 
size  distribution  of  the  channel  bed  sediment. 

4. 2. 3 . 4  At-A-Station  Downstream  Hydraulic  Geometry 

In  section  1.1  the  difference  between  downstream  and 
at-a-station  hydraulic  geometry  was  explained  as  a 
difference  arising  from  the  consideration  of  the  downstream 
variation  in  the  channel  characteristics  individually  in  the 
former  case  as  opposed  to  considering  the  variation  in  the 
response  of  the  characteristics  cumulatively  in  the  latter. 
At-a-station  downstream  hydraulic  geometry  is  simply  the 
comparison  of  various  stations  along  the  length  of  a  river. 
The  object  of  this  section  is  not  to  offer  an  additional 
explanation  for  the  variation  in  channel  behaviour  down 
Beauty  Creek  Flats  but  rather  by  considering  each  of  the 
three  reaches  as  stations  (discussed  previously  on  p.  70)  to 
distinguish  graphically  the  channel  behaviour  in  each  reach 
through  the  use  of  the  b-f-m  diagram  of  Rhodes  (1977).  This 
diagram  is  a  triangular  plot  of  the  values  of  the  b,  f,  and 
m  exponents  and  was  developed  by  Rhodes  to  aid  in  the 
interpretation  of  the  variation  displayed  by  hydraulic 
geometry  relationships.  Rhodes's  b-f-m  diagram  is  presented 
below  as  figure  4-6,  p.  101.  The  interpretative  potential  of 
the  diagram  arises  from  Rhodes's  division  of  it  into  ten 
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Figure  4-6.  The  b-f-m  diagram  of  Rhodes  (1977,  fig.  7).  Rhodes’s 

divison  of  the  diagram  into  ten  channel  types  aids  the 
downstream  at-a-station  interpretation  of  channel  change. 

The  variation  in  channel  response  to  an  increasing  discharge 
down  Beauty  Creek  Flats  is  indicated  by  the  plotting 
position  of  the  upper,  middle  and  lower  reach.  Also  shown 
is  the  plotting  position  for  Langbein’s  theoretical  b,  f 
and  m  values  for  noncohesive  channels,  the  point  representing 
the  Sunwapta  River  at  Beauty  Creek  Flats  and  the  point 
representing  the  White  River  of  Fahnestock,  1963. 


channel  types  (numbered)  using  five  criteria  drawn  mainly 
from  Leopold  and  Maddock's  1953  publication.  These  criteria 
are  represented  by  the  five  lines  on  the  diagram.  The 
interpretative  implications  of  each  of  these  divisions  is 
discussed  in  detail  in  Rhodes's  publication. 

The  presentation  of  this  diagram  serves  two  purposes. 
Firstly,  as  mentioned  above  it  allows  a  graphical 
demonstration  of  the  variation  in  channel  behaviour  between 
the  three  reaches  of  Beauty  Creek  Flats.  Secondly,  by 
plotting  the  position  of  the  three  reaches  on  it  and 
comparing  the  channel  behaviour  predicted  by  the  diagram 
(see  Rhodes's  table  2,  p.  82)  to  the  actual  channel 
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behaviour  the  interpretative  power  of  the  diagram  can  be 
tested. 

Use  of  the  diagram  requires  that  the  b,  f  and  m  values 
sum  to  unity  as  dictated  by  the  identity  Q  =  wdv.  However, 
in  the  majority  of  cases  they  do  not  (see  Ehodes,  1977,  p. 
75) .  In  order  to  plot  on  the  diagram  the  values  of  b,  f  and 
m  must  be  adjusted  so  that  they  do  equal  unity.  The 
adjustment  procedure  used  here  is  that  used  by  Ehodes  on  p. 
75.  Ehodes  felt  that  errors  in  measurement  of  depth  and  in 
particular  velocity  were  more  likely  than  errors  in  width 
measurement  and  therefore  adjusted  only  the  f  and  m 
exponents.  The  sum  of  the  b,  f  and  m  exponents  for  the 
upper,  middle  and  lower  reach  of  Beauty  Creek  Flats  are 
0.90,  0-92  and  0.96  respectively.  For  the  purpose  of 
plotting  on  the  b-f-m  diagram  the  f  and  m  values  are 
increased  by  0.05  each  in  the  upper  reach,  f  is  incremented 
by  +0.03  and  m  by  +0.05  in  the  middle  reach  and  both  f  and  m 
are  increased  by  0.02  in  the  lower  reach.  The  change  in  the 
plotting  position  of  the  adjusted  exponents  relative  to  the 
unadjusted  exponents  is  very  small  but  as  Ehodes  (1977,  p. 
75)  emphasized,  since  adjustments  were  made  the  positions  of 
the  points  can  only  be  considered  as  close  approximations  of 
the  correct  locations.  Included  on  the  diagram  is  the 
plotting  position  for  noncohesive  channels  determined  by 
using  the  theoretically  derived  b  (0.48),  f  (0.30)  and  m 
(0.22)  values  of  Langbein  (1965,  p.  306),  the  plotting 
position  of  the  White  Eiver  channels  (Fahnestock,  1963)  and 
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the  plotting  position  for  Beauty  Creek  Flats  whose  br  f  and 
m  values  are  0.42,  0.30  and  0.21  respectively.  For  plotting 
on  the  b-f-m  diagram  the  f  and  m  exponents  for  Beauty  Creek 
Flats  are  adjusted  by  +0.03  and  +0.04  respectively.  The 
White  River,  Beauty  Creek  Flats  and  Langbein's  theoretical 
b-f-m  values  plot  as  a  type  5  channel. 

Figure  4-6  shows  that  the  channels  of  the  upper  and 
middle  reach  plot  as  type  5  channels  while  the  channels  of 
the  lower  reach  plot  as  type  9.  As  such  the  points  lie  in 
the  region  of  the  diagram  which  Rhodes  considered  to  be 
indicative  of  braided  rivers  (see  Rhodes,  figure  5,  p.  86)  . 
Rhodes  (1977,  table  2,  p.  82)  has  tabulated  the  response 
predicted  for  a  type  5  and  type  9  channel  with  respect  to 
the  change  in  four  channel  parameters  with  increasing 
discharge.  Table  4-3  presented  below  compares  the  response 
of  these  parameters  to  an  increasing  discharge  as  predicted 
by  the  b-f-m  diagram  for  channel  types  5  and  9,  to  their 
actual  response  in  the  upper,  middle  and  lower  reach.  The 
bracketed  number  after  the  response  is  the  exponent  value 
determined  for  that  particular  power  function  relation.  As 
the  table  shows  the  predicted  response  of  the  channel 
parameters  is  exactly  that  which  actually  occurred  for  all 
three  reaches.  This  is  considered  to  be  strong  evidence  of 
the  interpretative  ability  of  the  b-f-m  diagram  in  hydraulic 
geometry  studies  of  channel  change. 

In  summary,  an  at-a-station  downstream  hydraulic 
geometry  plot  of  the  three  reaches  of  Beauty  Creek  Flats 
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indicates  a  change  from  a  type  5  channel  to  a  type  9  channel 
as  defined  by  Rhodes's  b-f-m  diagram.  The  plot  (see  figure 
4-6)  indicates  that  the  middle  reach  channels  although 
belonging  to  the  same  channel  type  as  the  channels  of  the 
upper  reach  are  distinguished  by  a  larger  b/f  ratio  and  a 
smaller  m  value.  The  lower  reach  channels  have  approximately 
the  same  b/f  ratio  as  the  upper  reach  and  are  distinguished 
primarily  by  a  low  rate  of  increase  in  velocity  with 
discharge,  or  in  other  words  a  significantly  smaller  m 
value. 

4. 2. 4  Best  Fit  Equations  For  Beauty  Creek  Flats 

Equations  4.13  to  4.16,  p.  76,  determined  for  the 
collated  Beauty  Creek  Flats  data  were  not  necessarily  the 
most  appropriate  mathematical  model  of  the  relationships 
between  the  dependent  variables  and  discharge.  As  stated  in 
section  1.1,  power  functions  are  the  conventional  model  used 
in  the  hydraulic  geometry  approach  to  channel  change.  To 
determine  if  other  models  would  provide  a  better  comparison 
between  the  observed  and  predicted  values  of  the  dependent 
variable,  polynomials  up  to  and  including  the  second  order 
were  tested.  Strictly  speaking  these  may  not  be  considered 
the  best  fit  equations  since  successively  higher  order  terms 
up  to  N-1  could  be  added  to  improve  the  predictive  ability 
of  the  model.  The  cutoff  at  the  second  order  term  was  based 
on  practicality  as  well  as  on  conceptual  grounds.  In  the 
latter  context,  since  the  dominant  discharge  was 
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approximately  bankfull  and  the  range  of  stage  during  the 

hours  of  the  day  when  the  data  was  gathered  was  reasonably 

limited,  not  even  quadratic  behaviour  would  be  expected  from 

any  of  the  dependent  variables-  This  expectation  was  proven 

to  be  correct  by  the  model  testing.  The  best  mathematical 

model  for  width,  average  depth  and  average  surface  velocity 

turned  out  to  be  the  power  functions  presented  as  equations 

4-13  to  4.15  inclusive,  p.  78,  while  that  for  water  surface 

slope  was  of  the  general  form  y  =  y.  That  this  is  the  best 

model  for  water  surface  slope  is  not  surprising  in  light  of 

2 

the  very  low  slope  (-0-0  1)  and  r  value  of  its  power 

function  (see  equation  4.16). 

Of  the  best  fit  models  the  most  statistically  sound  are 

equations  4.13  and  4.14.  The  choice  of  the  most 

statistically  sound  was  based  on  the  value  of  the 

2 

coefficient  of  determination  (r  )  and  on  the  appearance  of 

the  residuals  plot.  To  be  accepted  the  residuals  plot  had  to 

2 

display  a  roughly  tabular  pattern  and  the  r  values  had  to 
exceed  0-80  at  a  5%  level  of  significance.  Only  the  width 
(eg.  4.13)  and  average  depth  (eg.  4.14)  equations  met  these 
criteria.  The  standard  error  for  these  equations,  in  log 
units,  is  0.108  and  0.  069  respecti vely. 

Since  the  models  chosen  are  both  power  functions  it 
would  seem  worthwile  to  briefly  examine  the  conceptual  and 
technical  difficulties  associated  with  their  usage  and  in 
particular  to  determine  the  applicability  of  these 
difficulties  to  the  gravel  channels  examined  on  Beauty  Creek 
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Flats.  These  difficulties  were  discussed  by  Thornes  (1977). 
They  are  (1)  indeterminacy  of  channel  response  (2) 
nonequilibrium  flow  conditions  (3)  discontinuity  in  response 
and  (4)  nonlinearity  in  response.  The  former  two  are  not 
unique  to  power  functions  but  rather  are  conceptual 
difficulties  associated  with  channel  change  and  river 
hydraulics  in  general.  One  must  be  aware  of  them  but  their 
solution  awaits  the  development  of  a  more  precise  and 
thorough  theoretical  understanding  of  open  channel  flow.  For 
the  latter  two,  the  problem  associated  with  the  use  of  power 
functions  arises  from  their  log-linear  representation  of  the 
river  channel  relations-  Usually  this  is  an 
oversimplification  of  the  relationships  as  discontinuities 
or  nonlinear  response  in  the  dependent  parameter  is  often 
observed.  With  reference  to  discontinuity,  it  was  explained 
that  a  break  in  the  roughness-discharge  plots  (a  similar 
break  would  therefore  be  expected  for  velocity  and  depth 
-discharge  plots)  occurs  at  the  change  from  subcritical  to 
supercritical  flow  which  is  attributable  to  a  sudden 
reduction  and  then  increase  in  form  resistance  associated 
with  the  destruction  and  creation  of  dune  bed  forms.  Such 
scour  activity  is  usually  restricted  to  sand  bed  channels 
and  as  such  discontinuity  in  response  is  considered  to  be 
nonapplicable  to  the  gravel  channels  of  Beauty  Creek  Flats. 
Most  hydraulic  geometry  data  is  gathered  from  subcritical 
flows  and  therefore  such  discontinuities  usually  are  not 
important  considerations.  Nonlinearity  in  response  results 
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from  changes  in  the  rate  of  increase  or  decrease  (slope)  of 
the  dependent  variable  due  to  either  a  sudden,  drastic 
variation  in  flow  conditions  or  abrupt  changes  in  the  cross 
sectional  geometry  of  the  channel  with  rising  or  falling 
stage  (see  Richards,  1976)  or  both-  The  nonapplicability  of 
this  situation  to  the  collated  Beauty  Creek  Flats  data  has 
already  been  demonstrated  by  virtue  of  the  fact  that  second 
order  polynomials  did  not  fit  this  data  significantly  better 
than  the  log-linear  power  functions.  As  explained  above, 
this  is  the  result  of  the  lack  of  very  large  fluctuations  in 
stage  during  the  time  of  the  day  that  data  was  acguired;  ie- 
fluctuations  which  would  approximate  one  half  or  more  of 
average  channel  depth. 

In  summary,  as  a  result  of  the  coarse  nature  of  the 
Beauty  Creek  Flats  channels  and  due  to  the  dominant 
discharge  range  (approximately  bankfull)  at  the  time  of  data 
acquisition  discontinuities  or  nonlinearity  in  the  response 
of  the  dependent  parameters  did  not  occur.  The  result  is 
that  the  log-linear  power  function  representation  of  the 
relationships  do  not  constitute  severe  oversimplifications. 
Therefore,  their  utility  as  estimators  of  channel  response 
to  varying  discharge  is  at  a  maximum  under  these  flow 


conditions. 
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4-3  Comparison  Of  The  Beauty  Creek  Flats  Channels  To 
Channels  From  Other  Gravelly  Glacio fluvial  Environments 

4.3.  1  Introduction 

In  classifying  hydraulic  geometry  data  sets  according 
to  climate  Park  (1977)  included  the  studies  of  Nordseth 
(1973)  and  Brush  (1961)  in  the  proglacial  category.  However, 
reference  to  these  publications  indicates  that  the  rivers  of 
neither  are  currently  proglacial  and  are  therefore  not 
comparable  to  the  Sunwapta  River.  The  most  recent  major 
study  of  fluvial  processes  in  the  proglacial  environment  was 
done  by  Church  (1972)  on  Baffin  Island  sandar.  The  mean 
hydraulic  geometry  exponents  for  width,  depth  and  average 
velocity  on  these  sandar  are  provided  in  table  25,  p.  186  of 
Church's  publication.  Also  included  in  this  table  were  the 
width,  depth  and  average  velocity  exponents  from  Arnborg's 
1955  study  of  the  Hof f ellssandur  (computed  by  Church  from 
Arnborg's  data)  and  the  same  exponents  from  Fahnestock's 
1963  study  of  the  White  River,  Mount  Rainier,  Washington, 
U.S.A.  Although  they  are  all  studies  of  the  gravelly, 
glaciof luvial  environment  the  former  two  differ  from  the 
latter  one  with  respect  to  the  manner  in  which  the  channel 
data  was  acquired.  Church  and  Arnborg's  exponents  were 
determined  from  successive  flow  measurements  at  one 
(Arnborg)  or  a  number  of  gauging  stations  on  the  sandur(s) 
(Church)  and  as  such  represent  at-a-station  hydraulic 
geometry  exponents.  In  contrast  to  this  the  hydraulic 
geometry  exponents  of  Fahnestock's  study  were  derived  from 
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single  flow  measurements  at  many  different  stations  on  the 
White  River  valley  train.  As  was  mentioned  previously, 
Fahnestock  (1963,  p.  13)  indicated  that  this  approach  could 
be  considered  analogous  to  gathering  channel  data  from  a 
number  of  cross  sections  in  the  vicinity  of  a  particular 
gauging  station.  As  such  this  methodology  may  be  considered 
as  a  cruder  form  of  the  at-a-station  hydraulic  geometry 
approach  to  channel  change. 

The  methodology  used  to  acquire  the  channel  data  for 
the  Sunwapta  Eiver  at  Beauty  Creek  Flats  was  the  same  as 
that  used  by  Fahnestock  and  for  this  reason  only  his  study 
is  strictly  comparable. 

4-_J._2  Comparison  To  The  Channels  Of  The  White  Eiver 

Fahnestock,  over  a  period  of  two  years  (1958  and  1959) 
collected  data  on  112  channels  from  an  approximately  1.3  km 
stretch  of  the  White  River  between  the  confluence  of  west 
and  east  Emmons  Creek  and  the  point  where  the  river  is 
restricted  to  a  single  channel  by  a  moraine.  In  addition  to 
being  smaller  than  Beauty  Creek  Flats  (9  km)  the  White  River 
study  area  did  not  have  any  major  tributaries  (see 
Fahnestock's  figure  1  and  3,  p.  3  and  7).  In  order  to  have 
as  thorough  a  comparison  as  possible  ten  parameters  from 
Fahnestock's  table  7,  p.  62  (grain  size  data  taken  from 
table  2,  p.  24)  were  converted  to  metric  and  reduced  to 
their  summary  statistics.  This  data  is  presented  below  as 
table  4-4,  p-  111.  Comparison  of  this  data  to  the  data  on 


parameter 

X 

std.  dev. 

min/max 

N 

2 

flow  area  (m  ) 

1.70 

1.36 

0.002/5.64 

112 

S 

w 

0.023 

0.014 

0.002/0.062 

69 

w  (m) 

5.62 

3.28 

0.21/18.29 

112 

d  (m) 

0.24 

0.14 

0.01/0.63 

112 

w:d 

24.76 

10.28 

10/71 

112 

CO 

CO 

X 

l> 

1.12 

0.59 

0.09/2.74 

112 

Fr 

0.72 

0.23 

0.24/1.51 

112 

n 

0.057 

0.080 

0.013/0.085 

69 

D50  (mm) 

158.38 

24.98 

120/180 

8 

Q  (m  /s) 

2.55 

2.51 

0.0003/12.18 

112 

Table  4-4.  Statistical  summary  of  the  White  River  channel 
characteristics  of  Fahnestock’s  1963  study. 


the  60  channels  of  Beauty  Creek  Flats  contained  in  table  4-1 
indicates  a  high  degree  of  similarity  between  the  two  areas. 
Of  the  ten  parameters  compared  the  mean  values  for  area, 
width,  width  to  depth  ratio  and  average  velocity  on  the 
Sunwapta  exceeded  the  corresponding  values  for  the  White 
River.  For  water  surface  slope,  Froude  number,  roughness, 
discharge  and  grain  size  (D^ )  the  average  value  of  the 
White  River  was  greater  than  its  counterpart  on  the 
Sunwapta.  The  greatest  difference  between  the  two  areas  lies 
in  the  mean  values  for  slope  (0.023  vs  0.008),  Froude  number 
(0.72  vs  0.66)  and  grain  size  (White  River  vs  Sunwapta 
Fiver) .  It  should  be  noted  that  the  average  velocity  data 
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for  the  White  River  is  the  average  cross  section  velocity 
acquired  during  discharge  measurements  whereas  the  average 
velocity  for  the  Sunwapta  is  a  surface  velocity.  Also,  a 
strictly  valid  comparison  is  not  possible  for  grain  size 
since  his  D data  (see  his  table  2,  p.  24)  on  the  study 
area  mentioned  above  consisted  of  only  eight  data  points 
obtained  by  pebble  counts  of  material  greater  than  4  mm  in 
diameter  (the  method  of  Nolman,  1954,  was  used)  located  on 
subaerial,  central  portions  of  the  valley  train.  As 
mentioned  previously,  the  grain  size  data  for  Eeauty  Creek 
Flats  represents  grab  samples  taken  from  the  channel  bed. 

The  work  of  Ballantyne  (1975,  1978)  on  an  Ellesmere  Island 
sandur  and  Smith  (1974)  on  the  valley  train  of  the  Kicking 
Horse  River  in  Eritish  Columbia  indicates  that  the  channel 
deposits  are  coarser  than  adjacent  bar  sediment.  This 
suggests  that  despite  the  difference  in  sampling  technigue 
and  sample  location  the  D  ^  of  the  White  River  channels  are 
considerably  coarser.  Also  supporting  this  inference  are 
Fahnestock's  comments  on  boulder  movement  in  the  remarks 
column  of  his  table  7.  No  such  movement  or  sediment  occurred 
on  Beauty  Creek  Flats. 

The  greatest  similarity  between  the  two  rivers  is  found 

in  the  mean  values  for  the  average  depth  (0.24  vs  0.23  m) 

3 

and  discharge  data  (2.55  vs  2.47  m  /s;  White  vs  Sunwapta). 

An  hypothesis  test  for  the  comparison  of  two  sample 
means  using  the  discharge  means  of  the  two  areas  resulted  in 
the  acceptance  of  the  null  hypothesis  (equal  population 
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means)  at  an  alpha  level  of  0.005.  This  indicates  that  the 
discharge  data  of  the  White  River  could  have  been  obtained 
from  Beauty  Creek  Flats  of  the  Sunwapta  River  and  vice 
versa.  They  are  statistically  the  same  population. 

In  addition  to  being  similar  with  respect  to  the  ten 
parameters  mentioned  above  the  two  areas  are  also  very 
comparable  in  terms  of  most  of  their  hydraulic  geometry 
exponents  (slope  exponents  excepted;  see  table  4-5).  Of  the 
five  relations  for  the  White  River,  three  (width,  depth  and 
velocity)  were  determined  by  Fahnestock  and  two  (roughness 
and  slope)  were  determined  from  his  data  in  table  7  for  the 
purpose  of  comparison. 

The  similarity  of  the  exponents  in  light  of  the  large 
difference  in  average  slopes  between  the  two  areas  suggests 
that  slope  is  not  an  important  determining  parameter  of 
channel  behaviour  in  these  two  areas.  This  also  suggests 
that  slope  may  have  been  an  independent  parameter  in  the 
White  River  channels  as  it  has  been  interpreted  to  be  in  the 
channels  of  Beauty  Creek  Flats. 


4.3.3  Comparison  To  Channel  Behaviour  In  Other  Gravelly 
Glaciofluvial  Areas 

Table  4-5  presented  below  gives  the  hydraulic  geometry 
exponents  determined  for  the  Sunwapta  River  at  Beauty  Creek 
Flats  and  those  determined  for  the  White  River  (Fahnestock, 
1963) ,  Baffin  Island  Sandar  (Church,  1972)  and  the 
Hoffellssandur  (Arnborg,  1955).  Examination  of  this  table 
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exponent 

source 

W 

d 

V 

n 

s 

w 

Baffin  Island  sandar. 
Church  (1972) 

0.22 

0.31 

0.48 

/ 

/ 

Hoffellssandur,  Arnborg 
(1955) 

0.17 

0.39 

0.44 

/ 

/ 

White  River,  Fahnestock 
(1963) 

0.38 

0.33 

0.27 

0.01 

0.13 

Sunwapta  River  at 

Beauty  Creek  Flats 

0.42 

0.30 

0.21 

-0.01 

-0.01 

Table  4-5.  Comparison  of  the  exponents  of  the  hydraulic  geometry 

relations  determined  for  the  Sunwapta  River  at  Beauty  Creek 
Flats  to  the  corresponding  exponents  of  the  hydraulic  geometry 
relations  determined  for  similar  glaciofluvial  enviroments . 


indicates  that  for  the  gravel  channels  of  the  Baffin  Island 
Sandar  and  the  Hoffellssandur  the  primary  adjustment  to  an 
increasing  discharge  is  a  rapid  increase  in  velocity, 
whereas  for  the  gravel  channels  of  this  study  and  those  of 
the  White  River  the  primary  adjustment  to  an  increasing 
discharge  is  an  increase  in  channel  width.  Reference  to 
Rhodes's  (1977)  figure  5,  p.  86  indicates  that  the  channel 
behaviour  of  the  former  two  is  anomalous  for  braided 
channels  which  in  almost  all  cases  plotted  on  the  left  half 
of  the  b-f-m  diagram  (as  did  the  channels  of  the  Sunwapta 
and  White  River) .  This  is  a  significant  variation  in  channel 
behaviour  within  the  same  fluvial  environment.  Possible 


■' 
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reasons  for  this  variation  are  discussed  in  the  remainder  of 
this  section. 

In  order  to  be  able  to  increase  their  velocity  so 
rapidly  the  flow  resistance  in  the  channels  of  Church  and 
Arnborg's  studies  must  decrease  very  sharply  with  an 
increasing  discharge  (Church  determined  that  the  exponent  of 
the  Darcy-Weisbach  friction  factor-discharge  power  relation 
was  equal  to  -0.65).  Church  attributes  this  rapid  decrease 
primarily  to  a  rapid  increase  in  the  sediment  load  at  high 
discharges  which  would  smooth  out  the  roughness  of  the 
channel  bed  and  reduce  turbulence  through  an  increase  in  the 
viscosity  of  the  sediment-water  mixture.  Unfortunately,  the 
rate  of  increase  of  sediment  load  was  not  investigated  in 
this  study  and  is  not  known  for  either  Fahnestock's  or 
Arnborg's  study.  As  such  it  cannot  be  determined  whether 
this  can  account  for  the  variation  in  channel  behaviour 
between  the  two  groups.  However,  one  would  expect  that  in 
any  gravelly  proglacial  braided  river  a  relatively  rapid 
increase  in  sediment  load  would  accompany  the  rapid  increase 
in  discharge  that  characterizes  these  rivers.  It  is 
suggested  that  some  other  factor  which  is  common  to  Church's 
and  Arnborg's  studies  but  not  to  this  and  Fahnestock's  is 
responsible  for  the  variation  in  channel  behaviour.  This 
common  factor  is  suggested  to  be  the  low  temperatures  of 
high  arctic  summer  nights.  Church's  figure  15  indicates  that 
freezing  temperatures  frequently  occurred  during  June,  July 
and  August  at  Ekalugad  Fiord  and  from  data  gathered  at  three 
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climatological  stations  on  east  central  Baffin  Island  from 
1959  -  1968  he  determined  that  June's  mean  temperature  may 
have  exceeded  zero  degrees  Celsius  only  approximately  25  % 
of  that  time  (p.  13).  Table  4-5  indicates  that  the 

difference  between  the  Arctic  and  non-Arctic  studies  is  the 
magnitude  of  the  width  and  velocity  exponents,  the  depth 
exponent  being  similar  in  all  four  studies.  The  Arctic 
channels  experienced  much  lower  rates  of  change  of  width  and 
much  higher  rates  of  change  of  velocity-  If  freezing 
nocturnal  temperatures  are  to  account  for  this  anomalous 
behaviour  of  the  Arctic  channels  it  must  be  able  to  explain 
the  changes  in  these  exponents.  Low  width  and  high  velocity 
exponents  indicate  a  higher  bank  stability  which  the 
freezing  temperatures  must  therefore  be  able  to  explain.  Why 
freezing  temperatures  would  not  affect  the  depth  exponent 
must  also  be  explained. 

It  is  suggested  that  a  greater  bank  stability  is  caused 
by  freezing  of  intergranular  films  and  pockets 
(interstitial)  of  water  in  the  exposed,  wet,  bank  sedimer.t 
during  the  night  when  the  stage  has  dropped.  Thawing  during 
the  following  day  would  often  be  incomplete  due  to  a  very 
small  temperature  variation  (also  indicated  by  Church's 
figure  15  and  mentioned  on  his  p.  13) .  Incomplete  thawing 
would  be  aided  by  the  low  thermal  conductance  of  the  gravel. 
The  result  would  be  banks  which  are  more  resistant  to 
erosion  than  non-cohesive  banks  in  southern  latitudes.  This 
would  cause  a  low  rate  of  increase  of  channel  width  and 


channels  unable  to  accomodate  an  increase  in  discharge  by 
increasing  their  width  will  either  increase  their  depth  or 
velocity  in  an  attempt  to  maintain  within  channel  flow.  Of 
the  two,  an  increase  in  velocity  represents  the  least  work 
for  a  channel  flowing  in  gravelly  alluvium  resulting  in  it 
being  the  primary  adjustment  to  an  increased  discharge  in 
these  Arctic  channels.  The  channel  banks  of  this  and 
Fahnestock's  study  were  not  subjected  to  this  freezing  and 
as  a  result  their  primary  adjustment  is  an  increase  in  width 
due  to  the  ease  of  gravel  channel  bank  collapse. 

The  rate  of  change  of  depth  of  the  Arctic  channels  is 
not  affected  by  the  freezing  temperatures  due  to  the 
insulating  effect  of  the  channel  water  which  prevents  the 
channel  bed  sediment  from  freezing.  This  accounts  for  the 
similarity  in  the  depth  exponents  among  the  four  studies  of 
table  4-5. 

The  similarity  of  the  hydraulic  geometry  exponents 
between  this  study  and  Fahnestock's  and  between  Church's  and 
Arnborg's  studies  suggests  that  gravelly  proglacial  braided 
channels  behave  similarly  provided  the  associated  climatic 
conditions  are  similar-  It  is  therefore  reasonable  to 
suggest  equations  for  estimative  usage  on  areas  of  gravelly 
braided  outwash  other  than  those  included  in  table  4-5.  Such 
equations  could  also  be  used  for  estimating  the  paleoflow 
character istics  of  the  channels  in/on  exposed,  ancient, 
gravelly  braided  outwash. 

Since  Church's  study  is  the  only  major  study  of  braided 


' 
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gravel  outwash  in  the  Arctic  climate  his  width  and  velocity 
equations  are  suggested  for  estimating  channel  response  to 
an  increasing  discharge  in  this  climate.  In  the  non-Arctic 
climate,  the  only  two  are  this  study  and  that  of  Fahnestock 
(1963).  The  choice  of  equations  for  estimative  usage  from 
these  two  studies  is  based  on  the  value  of  the  coefficient 
of  determination.  In  this  context,  the  width  and  velocity 
equations  of  Fahnestock  were  chosen. 

Since  the  depth  response  did  not  vary  significantly 
between  the  Arctic  and  non-Arctic  studies  the  depth  relation 
with  the  best  coefficient  of  determination  was  chosen  for 
estimative  usage  in  both  climates.  The  depth  relation  of 
this  study  was  chosen  on  this  basis. 


4.4  Conclusions 

The  investigation  of  the  active  gravel  channels  on 
Beauty  Creek  Flats  of  the  Sunwapta  River  has  resulted  in  the 
following  conclusions  which  have  been  grouped  into  three 
categories. 

Conclusions  Applicable  To  Gravelly  Braided  Outwash  In 
General 


1.  There  is  a  similarity  in  channel  response  to  an 

increasing  discharge  in  the  gravelly  glaciof luvial 
environment  provided  the  associated  climatic  conditions 
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(specifically  temperature)  are  similar- 

2.  It  is  suggested  that  the  larger  velocity  exponent  and 
lover  width  exponent  of  braided  gravel  outwash  in  an 
Arctic  climate  may  be  due  to  a  greater  bank  stability 
caused  by  freezing  of  the  channel  bank  sediment  during 
the  night  when  the  stage  has  dropped. 

3.  The  width  and  velocity  equations  of  Church  (1972,  p. 

186)  and  the  depth  relation  of  this  study  are  suggested 
for  estimating  gravelly  glaciof lu vial  channel  response 
in  an  Arctic  climate.  In  the  non-Arctic  climate  the 
depth  relation  of  this  study  (eg.  4.14)  and  the  width 
and  velocity  equation  of  Fahnestock  (1963)  are  suggested 
for  similar  usage. 

4.  The  rate  of  change  of  depth  with  an  increasing  discharge 
in  braided  gravel  outwash  is  similar  in  Arctic  and 
non-Arctic  climates.  The  depth  relation  of  this  study  is 
suggested  for  estimative  usage. 

5.  The  rate  of  increase  of  width  with  increasing  discharge 
(b  exponent)  may  be  used  as  a  braiding  intensity  index. 

Conclusions  From  Channel  Data  On  Beauty  Creek  Flats 

1.  A  variation  in  channel  behaviour  does  exist  between  an 
upper,  middle  and  lower  reach  of  Beauty  Creek  Flats. 


2. 


The  most  influential  parameter  affecting  this  variation 
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is  the  median  grain  size  (D )  of  the  channel  bed  and 
banks  which  decreases  downstream. 

3.  In  all  three  reaches  of  Beauty  Creek  Flats  the  primary, 
secondary  and  tertiary  adjustment  of  the  channels  to  an 
increase  in  discharge  is  an  increase  in  width,  depth  and 
velocity  respectively. 

4.  Bankfull  flow  exceeds  the  critical  condition  required 
for  significant  gravel  movement  in  all  three  reaches. 

5.  The  highest  rate  of  sediment  deposition  occurs  in  the 
channels  of  the  middle  reach. 

6.  The  alluvial  fan  at  the  lower  end  of  Beauty  Creek  Flats 
exerts  a  dominant  control  only  over  the  velocity- 
discharge  relation. 

7.  The  upstream  growth  of  the  grassy  area  located 

I 

immediately  below  the  lower  reach  has  been  caused  by 
more  freguent  overbank  activity  in  this  area  of  Beauty 
Creek  Fiats. 

8.  Under  the  conditions  of  this  channel  study  (gravel  bed 
and  banks  and  a  limited  fluctuation  in  stage  around  the 
bankfull  condition)  power  functions  are  good  predictors 
of  channel  behaviour. 

9.  Rhodes's  (1977)  b-f-m  diagram  provided  an  accurate 
prediction  of  channel  response  to  an  increasing 
discharge  for  the  channels  of  this  study. 

10.  The  channels  of  the  middle  reach  are  distinguished  from 
the  upper  reach  channels  by  a  larger  b/f  ratio  and  a 
smaller  m  value.  The  lower  reach  channels  are 
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distinguished  from  those  of  the  upper  and  middle  reach 
by  a  significantly  smaller  m  value. 


Conclusions  From  Observations  On  Beauty  Creek  Flats 


1.  The  channels  in  the  active  flow  areas  of  Beauty  Creek- 
Flats  are  very  unstable.  The  major  cause  of  braiding  is 
erodible  banks. 

2-  Significant  gravel  movement  on  Beauty  Creek  Flats  occurs 
during  only  two  periods  of  the  year-  These  are  in  the 
spring,  during  runoff,  and  in  the  summer  during  peak- 
discharge. 

As  a  whole.  Beauty  Creek  Flats  is  neither  aggrading  or 
degrading.  Its  sediment  is  simply  being  reworked. 


3. 
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Appendix  A 

The  measured  discharge  data  set  (table  A-1),  the  theoretical 
discharge  data  set  (table  A-2)  and  table  A-3  are  presented 
in  this  appendix-  Tables  A-1  and  A-2  contain  the  data  used 
to  develop  the  hydraulic  geometry  relations  tor  the  Sunwapta 
River  at  Beauty  Creek  Flats.  The  details  of  the  field 
procedure  used  for  each  data  set  are  found  on  p.  43  -  47. 
With  reference  to  tables  A-1  and  A-2  several  points  should 
be  noted. 

1-  The  variability  of  the  gravel  to  sand  ratio  contained  in 
table  A-2  is  a  reflection  of  the  crudeness  of  the  shovel 
and  bucket  sampling  technique  employed.  However,  the 
decreasing  downriver  trend  of  this  ratio  shown  in  table 
4-1  demonstrates  that  the  channel  bed  samples  so 
obtained  are  still  representative  despite  the 
variability. 

2.  For  one  channel  (channel  6/1)  of  table  A-2  the  water 
surface  slope  was  zero  due  to  balancing  positive  and 
negative  gradients  along  the  channel.  In  this  situation 
the  bed  slope  was  used  for  any  calculations  requiring 
the  use  of  slope  data. 

3.  For  the  calculation  of  the  channel  Reynolds  number  for 
the  channels  of  the  measured  discharge  data  set  the 
average  kinematic  viscosity  for  the  upper  (transects  1 
and  2),  middle  (transects  3,  4  and  5)  and  lower  reach 
(transects  6  and  7;  see  figure  1-3)  of  the  theoretical 
discharge  data  set  was  determined  and  used. 
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4.  The  flow  area  of  table  A-1  represents  the  quotient  of 
discharge  over  the  average  flow  cross  section  velocity 
as  determined  by  current  meter  measurements.  For  the 
channels  of  table  A-2  no  such  velocity  was  obtained  and 
the  flow  area  was  determined  by  dividing  discharge  by 
the  average  velocity  at  0.6  of  the  depth. 

The  standard  equations  used  for  the  calculation  of  the 

channel  Reynolds  number  (  Re  , )  and  Frouae  number  (Fr)  of 

ch 

table  A-1  and  A-2  and  for  the  average  bottom  shear  stress 
(t  )  and  Mannings  n  of  table  A-2  are  presented  below. 


x  =  Y  R  S 

n  -h 


0 


Eg.  A.  1 


where 


R  =  d 


S  =  S 


w 


Eq.  A. 2 


where 


V  =  V 


0.6 


R  =  d 


v  =  y/p 
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Fr 


V 

% 

(gd)2 


where 


n 


r0.67  s0.50 
V 


E  q.  A  -  3 


Eq.  A.  4 


where  R  =  d 

S  =  S 

w 

V  =  V 

xs 


Note:  The  average  cross  section  velocity  of  the  theoretical 
discharge  data  set  was  derived  from  the  average,  average 
surface  velocity  to  average  cross  section  velocity  ratio 
(1.67)  determined  from  the  data  of  the  measured  discharge 


data  set. 


TABLE  A  -  1.  THE  MEASURED  DISCHARGE  DATA  SET 
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Appendix  B  -  Flow  Cross  Sections  Of  The  Theoretical 

Discharge  Data  Set 

Each  of  the  following  flow  cross  sections  has  an  identifying 
number  and  scale  presented  with  it-  The  cross  section  was 
identified  by  the  number  of  its  transect  (see  figure  1-3) , 
the  number  given  to  its  channel  on  that  transect  and  by  the 
position  of  the  cross  section  along  the  channel  length 
(upper,  middle  or  lower)-  For  example,  the  upper  flow  cross 
section  of  channel  number  two  on  transect  number  one  would 
be  given  the  identifying  number  1/2/T.  The  following  two 
points  should  be  kept  in  mind  when  examining  the  flow  cross 
sections : 

1.  All  the  flow  cross  sections  were  drawn  with  the 
assumption  of  a  horizontal  water  surface. 

2.  The  flow  cross  sections  were  constructed  with  the 
perspective  of  looking  south  from  highway  93- 
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Appendix  C  -  Grain  Size  Data  Of  The  Channel  Grab  Sanples 
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